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Abstract 

In an arbitrary unitary 4D CFT we consider a scalar operator 4>, and the operator 4> 2 
denned as the lowest dimension scalar which appears in the OPE <f> x <j) with a nonzero co- 
efficient. Using general considerations of OPE, conformal block decomposition, and crossing 
symmetry, we derive a theory-independent inequality [(f) 2 ] < f ([(/)]) for the dimensions of 
these two operators. The function f(d) entering this bound is computed numerically. For 
d — > 1 we have f(d) = 2 + 0(y/d — 1), which shows that the free theory limit is approached 
continuously. We perform some checks of our bound. We find that the bound is satisfied 
by all weakly coupled 4D conformal fixed points that we are able to construct. The Wilson- 
Fischer fixed points violate the bound by a constant 0(1) factor, which must be due to the 
subtleties of extrapolating to 4 — e dimensions. We use our method to derive an analogous 
bound in 2D, and check that the Minimal Models satisfy the bound, with the Ising model 
nearly-saturating it. Derivation of an analogous bound in 3D is currently not feasible be- 
cause the explicit conformal blocks are not known in odd dimensions. We also discuss the 
main phenomenological motivation for studying this set of questions: constructing models of 
dynamical Electro Weak Symmetry Breaking without flavor problems. 
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Figure 1: The best current bound (1.4), obtained by the method described in Section 5 



The subscript in / 6 refers to the order of derivatives used to compute this bound. 



1 The problem and the result 

Operator dimensions in unitary Conformal Field Theories (CFT) are subject to important con- 
straints known as unitarity bounds. In the simplest case of a scalar primary operator <ft, the 
unitarity bound states thatfj] 

d=[<p}>i, (l.i) 

d = 1 is free. (1.2) 

This classic result invites the following question: What happens if d = 1 +e? In particular, is there 
any sense in which the CFT (or at least its subsector not decoupled from 0) should be close to 
the free scalar theory if d is close to 1? For instance, do all operator dimensions in this subsector 
approach their free scalar theory values in the limit d — > 1? The standard proof of the unitarity 
bound P does not shed light on this question. 

In this paper we will show that such continuity indeed holds for the operator '0 2 ', by which 
we mean the lowest dimension scalar primary which appears in the OPE of <p with itself: 

0(x)0(O) ~ (x 2 )- d (l + C|x| Amin 2 (O) + ...), C^0. (1.3) 

In free theory A min = [eft 2 ] = 2, and we will show that A min — > 2 in any CFT as d — > 1. More 
precisely, we will show that in any 4D CFT 

A min < f{d), (1.4) 

where f(d) is a certain continuous function such that /(l) = 2. We will evaluate this function 
numerically; it is plotted in Fig. [I] for d near 1. 



We stress that bound (1.4) applies to the OPE 0x of an arbitrary scalar primary (f>. However, 
since the function f(d) is monotonically increasing, the bound is strongest for the scalar primary 
of minimal dimension. 



1 Unless explicitly noted otherwise, all statements of this paper refer to D = 4 spacetime dimensions. 
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Our analysis will use only the most general properties of CFT, such as unitarity, OPE, con- 
formal block decomposition, and crossing symmetry. The resulting bound (1.4) is thus model 
independent. In particular, it holds independently of the central charge, of the spectrum of what- 
ever other operators which may be present in the CFT, and of the coefficients with which these 
operators may appear in the OPE (1.3). Our analysis is non-perturbative and does not assume 
that the CFT may be continuously connected to the free theory. 

We do not know of any 4D CFT which comes close to saturating the bound of Fig. [TJ We are 
not claiming that this bound is the best possible, and in fact we do know that it can be somewhat 
improved using our method and investing some more computer time. In spite of not being the 
best possible, the curve of Fig. [T] is a valid bound, and represents a necessary condition which 
should be satisfied in any unitary CFT. 

The paper is organized as follows. In Section [2] we explain the phenomenological motivations 
behind this question, which are related to the naturalness problem of the electroweak scale. In 
Section [3] we review the necessary CFT techniques. In Section [4] we derive a sum rule for the 
contributions of all primary fields (with arbitrary spins and dimensions) appearing in the <fi x <p 
OPE. In Section [5] we explain how the sum rule is used to derive the bound (1.4). In Section [6] we 
check our bound against operator dimensions in various calculable CFTs in D = 4 and 4 — e. We 
also present and similarly check an analogous bound in D = 2. In Section [7] we discuss to what 
extent our result in its current form addresses the phenomenological problem from Section |2j In 
Section H] we conclude and indicate future research directions. 



2 Phenomenological motivation 

The phenomenological motivation of our study is given by one declination of the hierarchy problem, 
which was lucidly discussed in a paper by Luty and Okui [2] (see also [3]). This section is to a 
significant extent a review of the discussion in that paper. The bulk of the paper is logically 
independent of this section, and the reader who is mainly interested in the formal aspects of our 
result may skip to Section |3j 

The issue of mass hierarchies in field theory can be conveniently depicted from a CFT view- 
point. Indeed the basic statement that a given field theory contains two widely separated mass 
scales A IR <C A uv already implies that the energy dependence of physical quantities at A m <C 
E <C Auv is small, corresponding to approximate scale (and conformal) invariance. In the case of 
perturbative field theories the CFT which approximates the behaviour in the intermediate mass 
region is just a free one. For instance, in the case of non-SUSY GUT's, Air and Auv are re- 
spectively the Fermi and GUT scale, and the CFT which approximates behaviour at intermediate 
scales is just the free Standard Model. From the CFT viewpoint, the naturalness of the hierarchy 
A IR <C Auvj or equivalently its stability, depends on the dimensionality of the scalar operators 
describing the perturbations of the CFT Lagrangian around the fixed point. In the language of 
the RG group, naturalness depends on the relevance of the deformations at the fixed point. If 
the theory possesses a scalar operator Oa, with dimension A < 4, one generically expects^] UV 



'See concrete examples in Section 2.1 below 
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physics to generate a perturbation 

£pert — cA uv Oa. , (2-1) 

corresponding roughly to an IR scale 

Air = c^Auv • (2.2) 

Absence of tuning corresponds to the expectation that c be not much smaller than 0(1). If 4 — A 
is 0(1) (strongly relevant operator) a hierarchy between A IR and Auv can be maintained only by 
tuning c to be hierarchically smaller than one. This corresponds to an unnatural hierarchy. On 
the other hand when 4 — A is close to zero (weakly relevant operator) a mass hierarchy is obtained 
as soon as both 4 — A and c are just algebraically smalQ For instance for 4 — A = c = 0.1 
the mass hierarchy spans 10 orders of magnitude. Therefore for a weakly relevant operator a 
hierarchy is considered natural. The hierarchy between the confinement and UV scale in Yang- 
Mills theory is an example in this second class, albeit a limiting on^J The only exception to 
the above classification of naturalness concerns the case in which the strongly relevant operators 
transform under some global approximate symmetry. In that case it is natural to assume that the 
corresponding c's be small, even hierarchically small. The stability of the hierarchy depends then 
on the dimension As of the scalar singlet (under all global symmetries) of lowest dimension. If 
4 — As <C 1 the hierarchy is natural. 

According to the above discussion, in the SM the hierarchy between the weak scale and any 
possible UV scale is unnatural because of the presence of a scalar bilinear in the Higgs field 
H^H which is a total singlet with dimension ~ 2. On the other hand in supersymmetric ex- 
tensions of the SM, such scalar bilinears exist but their coefficient can be naturally chosen to 
be small. In a general supersymmetric model the weak scale is then naturally generated either 
by a marginally relevant deformation (dynamical supersymmetry breaking) or simply by adding 
strongly relevant supersymmetry breaking deformations with small coefficients (soft supersymme- 
try breaking). Technicolor models are instead similar to the case of YM theory: at the gaussian 
fixed point there are no gauge invariant scalars of dimension < 4. 

As far as the hierarchy is concerned these extensions are clearly preferable to the SM. However 
as far as flavor physics is concerned the SM has, over its extensions, an advantage which is also a 
simple consequence of operator dimensionality. In the SM the flavor violating operators of lowest 
dimensionality, the Yukawa interactions, have dimension = 4, 

C Y = y%Hq L u R + yt j rfq L d R + yt j rfL L e R (SM) , (2.3) 

and provide a very accurate description of flavor violating phenomenology. In particular, the 
common Yukawa origin of masses and mixing angles leads to a critically important suppression 
of Flavor Changing Neutral Currents (FCNC) and CP violation. This suppression is often called 
Natural Flavor Conservation or GIM mechanism [4j. Once the hierarchy v <C Auv is taken 
as a fact, no matter how unnatural, extra unwanted sources of flavor violation are automatically 
suppressed. In particular the leading effects are associated to 4-quark interactions, with dimension 

3 We stole this definition from ref. [3]. 

4 This is because the corresponding deformation, the glueball field G^ V G^ , is marginally relevant: its scaling 
dimension is 4 — ag 2 and becomes exactly 4 at the gaussian fixed point. 
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6, and are thus suppressed by v 2 /A^j V . The situation is not as good in supersymmetry, where in 
addition to the Yukawa interactions flavor is violated by operators of dimension 2 and 3 involving 
the sfermions. The comparison with technicolor brings us to discuss the motivation for our paper. 
In technicolor the Higgs field is a techni-fermion bilinear H = TT with dimension ~ 3. The 
SM fermions instead remain elementary, i.e. with dimension 3/2. The Yukawa interactions are 
therefore irrelevant operators of dimension 6, 

C Y = ^Hqq (TC), (2.4) 

and are associated to some new dynamics [5], the flavor dynamics, at a scale Ap, which plays the 
role of our Auv- Very much like in the SM, and as it is found in explicit models [SJ, we also expect 
unwanted 4-quark interactions 

-rrQiQjQkqi (2-5) 

A F 

suppressed by the same flavor scale. Unlike in the SM, in technicolor the Yukawa interactions 
are not the single most relevant interaction violating flavor. This leads to a tension. On one 
hand, in order to obtain the right quark masses, Ap should be rather low. On the other hand, 
the bound from FCNC requires A-p to be generically larger. For instance the top Yukawa implies 



Af < lOTeV. On the other hand the bound from FCNC on operators like Eq. (2.5) is rather 
strong. Assuming Cijki ~ 1, flavor mixing in the neutral kaon system puts a generic bound ranging 
from Af > 10 3 TeV, assuming CP conservation and left-left current structure, to Ap > 10 5 TeV, 
with CP violation with left-right current. Of course assuming that c^-jy have a nontrivial structure 
controlled by flavor breaking selection rules one could in principle obtain a realistic situation. It 
is however undeniable that the way the SM disposes of extra unwanted sources of flavor violation 
is more robust and thus preferable. The origin of the problem is the large dimension of the Higgs 
doublet field H. Models of walking technicolor (WTC) |6J partially alleviate it. In WTC, above the 
weak scale the theory is assumed to be near a non-trivial fixed point, where H = TT has a sizeable 
negative anomalous dimension. WTC is an extremely clever idea, but progress in its realization 
has been slowed down by the difficulty in dealing with strongly coupled gauge theories in 4D. 
Most of our understanding of WTC relies on gap equations, a truncation of the Schwinger- Dyson 
equations for the TT self-energy. Although gap equations do not represent a fully defendable 
approximation, they have produced some interesting results. In case of asymptotically free gauge 
theory they lead to the result that H = TT can have dimension 2 at the quasi-fixed point, but not 
lower [7] . In this case the Yukawa interactions would correspond to dimension 5 operators, which 



are more relevant than the unwanted dimension 6 operators in Eq. (2.5). However some tension 
still remains: the top Yukawa still requires a Flavor scale below the bound from the Kaon system, 
so that the absence of flavor violation, in our definition, is not robust. It is quite possible that 
the bound [H] > 2 obtained with the use of gap equation will not be true in general. Of course 
the closer [H] is to 1, the higher the flavor scale we can tolerate to reproduce fermion masses, and 



the more suppressed is the effect of Eq. (2.5). However if [H] gets too close to 1 we get back the 
SM and the hierarchy problem! More formally, a scalar field of dimension exactly 1 in CFT is a 
free field and the dimension of its composite H^H is trivially determined to be 2, that is strongly 
relevant. By continuity we therefore expect that the hierarchy problem strikes back at some point 
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as [H] approaches 1. However the interesting remark made by Luty and Okui [2 J is that, after all, 
we do not really need [H] extremely close to 1. For instance [H] = 1.3 would already be good, 
in which case the corresponding CFT is not weakly coupled and it could well be that [wH\ is 
significantly bigger than 2 [if] and maybe even close to 4. The motivation of our present work 
is precisely to find, from prime principles, what is the upper bound on Aj = [WH\ as d — [H] 
approaches 1. In simple words this may be phrased as the question: how fast do CFTs, or better 
a subsector of CFTs, become free as the dimension of a scalar approaches 1? 

In the following subsection we would like to make a more quantitative analysis of the tension 
between flavor and electroweak hierarchy in a scenario where the electroweak symmetry breaking 
sector sits near a fixed point between the EW scale A IR and some A uv at which, or below which, 
Flavor dynamics must take place. 



2.1 Quantitative analysis 

Let us normalize fields and couplings in the spirit of Naive Dimensional Analysis (NDA). The 
Lagrangian will thus be written as 

C=^Ffa\,M), (2.6) 

where fields, couplings and physical mass scales are indicated collectively. With this normalization, 
Green's functions in the coordinate representation have no factors of ir, the couplings A are loop 
counting parameters, and the mass scales M correspond to physical masses (as opposed to decay 
constants). 

Our hypothesis is that below some UV scale Auv the theory splits into the elementary SM 
without Higgs and a strongly coupled CFT which contains a (composite) Higgs scalar doublet 
operator H . These two sectors are coupled to each other via weak gauging and the Yukawa 
interactions. As a warmup exercise consider then the top Yukawa 

-^\ t HQ L t R + h.c. (2.7) 
and its lowest order correction to the CFT action 

A£= (j^j X 2 J d 4 xd i yH(xyH(y)Q L t R (x)i R Q L (y) 



1 



16tt 2 
1 



2 

\ 2 t I d A xd A y{H^H){x)\x-yf s ~ 2d \x-y'- { ' 



16n 2 
1 



d A x\ 2 t Al + Y 2d - As (H^H)(x) 

d 4 x[\Mvv)] 2 Atjv As (H j H) (2. 



167T 2 

where we used the x H OPBjand cut off the d A y integral at a UV distance A^y (ex.: there exist 



5 Notice that H^H is defined as the scalar SU(2) singlet operator of lowest dimension in the W x H OPE. Weak 
SU(2) is assumed to be a global symmetry of the CFT. The SU(2) invariance of the fermion propagators realizes 
the projection on the singlet. 
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new states with mass A uv ). The quantity Ai(Auv) = AfA^y 1 represents the dimensionless running 
coupling evaluated at the scale Auv- Given our NDA normalization, A^ is the loop counting 



parameter (no extra 7r's). For As < 4 Eq. (2.8) represents a relevant deformation of the CFT. 



Allowing for a fine tuning e t between our naive estimate of the top loop and the true result, the 
deformation lagrangian is 



1 v 

67T 

corresponding to a physical infrared scale 



Ca =— A t (A uv ) 2 64- v A ^^ , (2.9) 

107T Z 



Air = [A t (A uv ) 2 e t ]^sA uv = (c t )^sA vv , (2.10) 

where we made contact with our previous definition of the coefficient c. If the above were the 
dominant contribution, then a hierarchy would arise for 4 — As < 1 provided At(Auv) < 1 and/or 
a mild tuning e t < 1 exists. However, unlike the normal situation where the Higgs is weakly self- 
coupled and the top effects dominate, in our scenario the Higgs is strongly self-coupled. Therefore 
we expect a leading contribution to Ca to already be present in the CFT independently of the 
top: 

C A = ^A^H^H. (2.11) 

The presence of such an effect basically accounts for the fact that a CFT with a relevant defor- 
mation does not even flow to the fixed point unless the deformation parameter is tuned. This is 
in line with our initial discussion. So we shall work under the assumption that c is somewhat less 
than 1. 

We can describe the generation of the electroweak scale by writing the effective potential for the 
composite operator vacuum expectation value ((H^H)) = fx As . Compatibly with scale invariance 
there will also be a term 

^CFT = ^ (2.12) 

where a is a numerical coefficient that depends on the CFT and on the direction of the VEV in 
operator space. The full effective potential has then the form 

V* = ^ h A m A V S + A > (2-13) 

which is stationary at /i ~ Ajr. Here we put for example a = 1 and chose a negative sign for the 
scale breaking contribution. Notice that the vacuum dynamics picture that we just illustrated is 
analogous to the Randall-Sundrum model |8] with Goldberger-Wise radius stabilization [9] with 
the identification of with the position of the IR brane in conformal coordinates. This fact 
should not be surprising at all given the equivalence of that model to a deformed CFT [TO] . 

Our analysis of the top sector is however useful to discuss the two basic constraints on this 
scenario. 

The first constraint is the request that Auv be below the scale where the top Yukawa becomes 
strong, at which point the SM becomes strongly coupled to the CFT and our picture breaks 
down. The running top coupling in its standard normalization is yt(E) = AttX(E). Using the 
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known experimental result ^(Atr) 

yt(Auv) 



1 we thus have 
'Auv 



y t (A 



A 



IR 



d-1 



Perturbativity corresponds to yt(A) < An, that is 

d- 1 



4- Ac 



ln(l/c) < ln(4vr) ~ In 10 . 



(2.14) 



(2.15) 



Clearly, this bound is better satisfied the closer d is to 1. 

A second constraint is presented by the request of robust decoupling of unwanted flavor break- 
ing effects. Assuming any generic interaction among SM states is present at the scale Auv with 
strength comparable to A 2 (A), we can parameterize flavor violation by 



C 



fermion 



167T 2 



which by going to canonical normalization becomes 

•^fermion 

with 



A 2 



Ap = A 



IR 



A 



uv 



A 



IR 



2-d 



Air(c) 



(2.16) 



(2.17) 



(2.18) 



By taking Air = 1 TeV, the bound from FCNC can be parametrized as Ap > 10 F TeV. Making 
the conservative assumption that all quark families appear in Eq. (2.16), compatibility with the 
data requires Ap > 10 3 4- 10 4 TeVjjThus robust suppression of FCNC corresponds to 



F > 3 4 4. (robust) . 



(2.19) 



On the other hand if only the third family appears in Eq. (2.16), the mixing effects involving the 
lighter generations are generally suppressed by extra powers of the CKM angles. In that case the 
bound on A F is weaker, and F > 0.5 is basically enough. In the latter case the detailed structure 
of the Flavor theory matters. Notice that for conventional walking technicolor models, for which 
d > 2, we always have Ap < Air, so that even the weaker bound is somewhat problematic. We 
are however interested to see to what extent we can neglect this issue by focussing on the robust 
bound (2.19). From Eq. (2.18) we must have 

2-d 



4- Ac 



ln(l/c) > Fin 10. 



Eqs. fl2.15|2.20p together imply 
F ln(l/c) 



< 



< 



2-d (4 -A s ) In 10 d-1 



d<l + 



l + F 



(2.20) 



(2.21) 



6 We consider the limit Ap = 10 5 TeV really an overkill. 
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which for the robust bound (2.19) requires d < 1.25. At the same time the amount of tuning 
needed to generate the hierarchy between Air and the flavor scale Auv is 



10 (4-A s )F/(d-2) 



A reasonable request c > 0.1 then reads 



4- A s < 



(2.22) 



The ultimate goal of our study is thus to find a prime principle upper bound on As as a function of 



d. This bound should provide an extra important constraint which together with eqs. ( 2.21|2.22 ) 
may or may not be satisfied. Our main result (1.4) is a step towards this goal, although is 



not yet a complete solution. The point is that the lowest dimension scalar in the 0x0 OPE, 



whose dimension A appears in (1.4), is not necessarily a singlet. Nonetheless we think our result 
already represents some interesting piece of information. We postpone a detailed discussion of 
this connection until Section 



3 Necessary CFT techniques 

To make the paper self-contained, in this section we will review a few standard CFT concepts and 
results, concentrating on those which are crucial for understanding our result and its derivation. 
Our personal preferred list of CFT literature includes [UJ , [TJ] , [13] , [H] , [15] , [16] . We will mostly 
work in the Euclidean signature. 

3.1 Primary fields and unitarity bounds 

In perturbative field theories, classification of local operators is straightforward: we have a certain 
number of fundamental fields, from which the rest of the operators are obtained by applying 
derivatives and multiplication. In CFTs, a similar role is played by the primary fields. These local 
operators 0(x) are characterized by the fact that they are annihilated by the Special Conformal 
Transformation generator (at x = 0). Thus a primary field 0(x) transforms under the little 
group=the subgroup of conformal transformations leaving x = invariant (this includes Lorentz 
transformations M M „, dilatations D, and Special Conformal Transformations K^) as follows: 

[M^, (9(0)] = £^(9(0), (3.1) 
[D,O(0)]=iAO(0), 

Here we assume that O has well-defined quantum numbers: the scaling dimension A, and the 
Lorent^jspin (the matrices are the corresponding generators). 

Once all the primary operators are known, the rest of the field content is obtained by applying 
derivatives; the fields obtained in this way are called descendants. The multiplication operation 

7 Or Euclidean rotation, if one is working in the Euclidean. 
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used to generate composite operators in perturbative field theories has a CFT analogue in the 



concept of the OPE, which will be discussed in Section pO] below. To avoid any possible confusion, 
we add that this picture applies equally well also to conformal gauge theories, e.g. to Af = 4 
Super Yang-Mills, provided that only physical, gauge invariant fields are counted as operators of 
the theory. 



Knowing (3.1), one can determine the transformation rules at any other point x using the 



conformal algebra commutation relations [llj. In principle, one could also imagine representa- 



tions where acts as a nilpotent matrix (type lb in [TT]) rather than zero as in (3.1). However, 



as proven in [T], only representations of the form (3.1) occur in unitary CFTs. Moreover, uni- 
tarity implies important lower bounds on the operator dimensions. We are mostly interested in 
symmetric traceless fields (/x) = Hi ■ ■ ■ Hh which correspond to j = j tensors: 

— a aicei • • • u a. x 6ii ( P ■ 

This is traceless in any pair of /i indices because cr^cr^ oc e^e^. For such primaries the unitarity 
bound reads pQ: 

I = 0: A > 1, A = 1 only for a free scalar; 

I > 1: A > I + 2, A = / + 2 only for a conserved current. (3-3) 

Notice a relative jump of one unit when one passes from I = to / > 1. In particular, a conserved 
spin-1 current has A = 3, while the energy- momentum tensor has A = 4. The full list of unitarity 
bounds, which includes also fields with j ^ j, can be found in pQ. Recently [TF], some of these 
bounds were rederived in a very physically transparent way, by weakly coupling a free scalar theory 
to the CFT and studying the unitarity of the S-matrix generated by exchanges of CFT operators. 

3.2 Correlation functions 

As is well known, conformal symmetry fixes the coordinate dependence of 2- and 3-point functions 
of primary fields. For example, for scalar primaries we have: 

W*)M)> = 1 ~~~|2a7 , (3-4) 

\x — y\ 4 

((f>(x)4>(y)) = (0^0). (3.5) 



As it is customary, we normalize to have a unit coefficient in the RHS of (3.4). Correlators of 
two fields with unequal dimensions vanish by conformal symmetry. Even if several primaries of 
the same dimension exist, by properly choosing the basis we can make sure that the nondiagonal 



correlators (3.5) vanish. Notice that we are working with real fields, corresponding to hermitean 
operators in the Minkowski space description of the theory. 
The 3-point functions are also fixed by conformal symmetry: 

(0i(^i)02(^2)03(^3)) Au " 



Xl2 |A 1 +A 2 -A 3 | a;23 |A 2 +A3-A 1 | Xl3 | 



X12 = x\ — X2 etc. 
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The constants A123, which become unambiguously defined once we normalize the fields via the 2- 
point functions, are an important characteristic of CFT dynamics. These constants appear as the 
OPE coefficients (see below), and if they are all known, any n-point function can be reconstructed 
via the OPE. Thus in a sense finding these constants, together with the spectrum of operator 
dimensions, is equivalent to solving, or constructing, the theory. 

Also the correlator of two scalars and a spin I primary Oi^ is fixed up to a constant 



7 _ ^13 _ x 23 
A 1 2 2 

x 13 x 23 



The OPE coefficients A120 are real, once a real field basis ( 3.4[ ) (and similarly for higher spin 



primaries) is chosen. This reality condition follows from the reality of Minkowski-space correlators 
of hermitean operators at spacelike separation, see Appendix [A] for a more detailed discussion. 

When it comes to 4-point functions, conformal symmetry is no longer sufficient to fix the 
coordinate dependence completely. In the case of 4 scalar operators, the most general conformally- 
symmetric expression is 

/1 I \ Ai — A2 /1 I \ A3 — A4 / \ 

(01(^1)02(^2)03(^3)04(^4)) = I I r I I ] 1 I 1 iaTZaTT IA.4-A, ' 

(3.7) 



v |a: 14 | y Vl x i3l/ |a;i2| Al+A2 |a;34| A 3 +A4 ' 

where g(u,v) is an arbitrary function of the conformally-invariant cross-ratios: 

22 22 

X 12 X 34 x 14 x 23 in o\ 

U = k-, V = k-. (3.8) 

J, 13 X 24 x 13 x 24 

3.3 Operator Product Expansion 

A very powerful property of CFT is the Operator Product Expansion (OPE), which represents a 
product of two primary operators at finite separation as a sum of local primaries: 

0i(x)0 2 (O) = ^A 12O [C (M) (x)(9 (M) (0) + ...], (3.9) 
o 

( x ) = j^jAi+Aa-Ac Up • 

Here we wrote an OPE appropriate for a pair of scalars 0i x <fi 2 - 



The ... in (3.9) stands for an infinite number of terms, less singular in the x — > limit, 
involving the derivatives of the primary O^) (i.e. its descendants). The coordinate dependence 
of the coefficients of these descendants is in fact completely fixed by the conformal symmetry, so 
that A120 appears as an overall coefficient for the full contribution of O'^ and its descendants. 
We can write schematically: 

01 X 02 = ^ A120 y— O , 
O 
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denoting by O the contribution of the whole conformal family. 

For example, for a scalar operator O appearing in the OPE <p x <ft the first few subleading 
terms are (pj], p. 125, A = A ) 

mm - j^gpx [i + \ x A + ~^x^d v - l^-^-^d 2 + ■ ■ wo) (3.io) 

There are several ways to determine the precise form of the coefficients of the descendants. 



One, direct, way [12] is to demand that the RHS of (3.9) transform under the conformal algebra 
in the same way as the known transformation of the LHS. A second way [25] is to require that the 
full OPE, with the descendant contributions included, correctly sum up to reproduce the 3-point 



function (3.6) not only in the limit X\ — ► X2, but at finite separation as well. The last, seemingly 
the most efficient way, is via the so called shadow field formalism, which introduces conjugate 
auxiliary fields of dimension 4 — A and uses them to compute "amputated" 3-point functions, 
which turn out to be related to the OPE coefficient functions [H]. 

Using the OPE, any n-point function (0!(x)0 2 (O) . • •) can be reduced to a sum of (n— l)-point 
functions. Applying the OPE recursively, we can reduce any correlator to 3-point functions which 
are fixed by the symmetry. Of course, this procedure can be carried out in full only if we already 



know which operators appear in the OPE, and with which coefficients. Consistency of (3.9) and 



(3.6) in the limit X\ —> 22 requires that the same constant X120 appear in both equations. Thus 



the sum in (3.9) is taken over all primaries 0(a) which have non-zero correlators (3.6). It is not 
difficult to show that the correlator ((frifoO) vanishes if O has j 7^ jj^Jand thus such fields do not 
appear in the OPE of two scalars (see e.g. [20], p. 156). 

We stress that in CFT, the OPE is not an asymptotic expansion but is a bona fide con- 
vergent power-series expansion^} The region of expected convergence can be understood using 
the state-operator correspondence in the radial quantization of CFT (see [IE], Sections 2.8,2.9 
for a lucid discussion in 2D). In this picture, every state \ty) defined on a sphere of radius r 
around the origin can be expanded in a basis of states generated by local operator insertions at 
the origin acting on the vacuum: O(0)|0). For example, consider the Euclidean 4-point function 
(0i(x)0 2 (O)03(x 3 )04(x4)), and suppose that 

< \x\ < min(|a;3|, |x 4 |) , (3-H) 

so that there exists a sphere centered at the origin which contains 0, x but not 23, X4. Cutting the 
path integral along the sphere, we represent the 4-point function as a Hilbert-space product 

(0 1 (a;)02(O)03(x3)04(x4)) = (*s\<hMM*d\0) , = (0\Mx)M0) , 

where the radial quantization state \^s) lives on the sphere, and can be expanded in the basis of 



local operator insertions at x — 0. Thus we expect the OPE to converge if (3.11) is satisfied. To 



Fields with j ^ j correspond to antisymmetric tensors. The correlator \ 4>x(x)(j)2(—x)0 ^^(0) ) must vanish for 
this particular spacetime configuration, since we cannot construct an antisymmetric tensor out of x^ L . Any other 
configuration can be reduced to the previous one by a conformal transformation. 

9 In general, it will involve fractional powers depending on the dimensions of the entering fields. 
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quote [IS], the convergence of the OPE is just the usual convergence of a complete set in quantum 
mechanics. See also [I9],[20] for rigorous proofs of OPE convergence in 2D and 4D CFT, based 
on the same basic idea. 

The concept of OPE is also applicable in theories with broken scale invariance, e.g. in asymp- 
totically free perturbative field theories, such as QCD, which are well defined in the UV. These 
theories can be viewed as a CFT with a relevant deformation associated to a scale Aqcd- The 
question of OPE convergence in this case is more subtle. In perturbation theory, the OPE pro- 
vides an asymptotic expansion of correlation functions in the x — > limit |21| . and is unlikely to 
be convergent because of non-perturbative ambiguities associated with the renormalons and the 
choice of the operator basis. It has however been conjectured in [22] that full non-perturbative 
correlators should satisfy a convergent OPE also in theories with broken scale invariance. This 
presumably includes QCD, but so far has been proved, by direct inspection, only for free massive 
scalar (see [21] and Ref. [6] of [22]). 



3.4 Conformal blocks 



As we mentioned in Section 3.2[ conformal invariance implies that a scalar 4-point function must 
have the form (3.7), where g{u,v) is an arbitrary function of the cross-ratios. Further information 
about g(u,v) can be extracted using the OPE. Namely, if we apply the OPE to the LHS of 



(3.7) both in 12 and in 34 channel, we can represent the 4-point function as a sum over primary 



operators which appear in both OPEs: 



^20304) — 2^ ^120-^340 



CB 



O 




(3.12) 
(3.13) 



The nondiagonal terms do not contribute to this equation because the 2-point functions of non- 
identical primaries O ^ O' vanish, and so do 2-point functions of any two operators belonging to 
different conformal families. The functions CB^, which receive contributions from 2-point func- 
tions of the operator O and its descendants, are called conformal blocks. Conformal invariance of 
the OPE implies that the conformal blocks transform under the conformal group in the same way 
as (01020304)- Thus they can be written in the form of the RHS of (3.7), with an appropriate 



function go(u,v). In terms of these functions, ( |3.12 ) can be rewritten as 



g{u,v) 



2^ ^120-^340 9o(u,v). 



(3.14) 



o 



In general, functions go(u,v) depend on the spin I and dimension A of the operator O, as 
well as on the dimensions Aj = [0;]. Various power-series representations for these functions were 
known since the 70's, but it seems that simple closed-form expressions were obtained only recently 
by Dolan and Osborn [25], [26]. In what follows we will heavily use their result, in the particular 
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case when all Aj are equal. In this case it takes the form independent of A, 
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(-)* zz 
z - z 

k p {x) =^/ 2 2 Fi(/5/2,/5/2,/5;x), (3.15) 



go(u,v) = gAj(u,v) = — - [k A+ i(z)k A -i-2(z) - (z <-> z)] , 



where the variables z, z are related to u, v via 

m = zz, v — (1 — z)(l — I), (3.16) 

or, equivalent lyF] 



- (it - V + 1 ± ^JJu 



z, z = - [u — v + 1 ± \/(u — v + l) 2 — Au 



We will give a brief review of the derivation of Eq. (3.15) in Appendix |Bj A short comment is 
here in order about the meaning and range of z and z (see Appendix [C] for a more detailed discus- 
sion). With points Xi varying in the 4D Euclidean space, these variables are complex conjugates 
of each other: z = z*. Configurations corresponding to real z = z can be characterized as having 
all 4 points lie on a planar circle. Below we will find it convenient to analytically continue to 
the Minkowski signature, where z and z can be treated as independent real variables. One possi- 
ble spacetime configuration which realizes this situation (the others being related by a conformal 
transformation) is shown in Fig. [2j Here we put 3 points along a line in the T = Euclidean 
section: 

Xl = (0,0,0,0), x 3 = (1,0,0,0), x 4 = oo, (3.17) 
while the 4-th point has been analitically continued to the Minkowski space: 

x 2 -> x ? = (X u 0,0,T), T = -iX A . (3.18) 

One shows (see Appendix [C]) that in this situation 

z = X 1 -T, z = X 1 + T. (3.19) 

We do not expect any singularities if x^ 1 stays inside the "spacelike diamond" region 

0<z,z<l, (3.20) 

formed by the boundaries of the past and future lightcones of X\ and x 3 . Indeed, one can check 
that the conformal blocks are real smooth functions in the spacelike diamond. 
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This happens because for equal Aj the coefficients of conformal descendants in the OPE are determined only 
by A and I, see e.g. Eq. fl3.10|). 



11 Notice that the RHS of (3.15 ) is invariant under 
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Figure 2: The "spacelike diamond" (3.20) in which the conformal blocks are real and 
regular, see the text. 



4 Crossing symmetry and the sum rule 

Let us consider the 4-point function ( |3.7[ ) with all 4 operators identical: 0j = 0. We have: 

g(u,v) 



\xi2\ 2d \x 3 A 2d ' 



d 



(4.1) 



The LHS of this equation is invariant under the interchange of any two Xi, and so the RHS should 
also be invariant, which gives a set of crossing symmetry constraints for the function g(u,v). 
Invariance under X\ <-> x 2 and X\ «-> £3 (other permutations do not give additional information) 
implies: 



g{u,v) = g(u/v, l/v) 
v d g{u,v) = u d g{v,u) 



(xi <-> x 2 ) 
(xi x 3 ) 



(4.2) 
(4.3) 



At the same time, g(u, v) can be expressed via the conformal block decomposition (3.14), which 
in the considered case takes the form: 



g(u,v) = l+ ^2 X o9o(u,v). 

0£<l>x<f> 



(4.4) 



Here in the first term we explicitly separated the contribution of the unit operator, present in the 



0x0 OPE. Since \q are real (see Section 3.2), all conformal blocks appear in (4.4) with positive 
coefficients. 



Let us now see under which conditions Eq. (4.4) is consistent with the crossing symmetry. The 



X\ <-> £2 invariance turns out to be rather trivial. Transformation properties of any conformal 
block under this crossing depend only on its spin [25J: 



9A,l[ u : v ) 



l g&,i{u/v,l/v). 
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All the operators appearing in the OPE 0x0 have even spir Thus the first crossing constraint 



(4.2) will be automatically satisfied for arbitrary coefficients \ 2 a . 



On the other hand, we do get a nontrivial condition by imposing that (4.4) satisfy the second 



crossing symmetry (4.3). This condition can be conveniently written in the form of the following 
sum rule: 



1 = ^2 P&J F dAJ ( z >z), VA,i > , 

A,l 

v d gA,i(u,v) - u d g A)l (v,u) 



Fd,A,l( z i z ) = 



(4.5) 



where the sum is taken over all A, I corresponding to the operators O G x 0, p^i = X%, and 



u, v are expressed via z,z via (3.16). As we will see below, this sum rule contains a great deal 



of information. It will play a crucial role in the derivation of our bound on the scalar operator 
dimensions. 



Below we will always apply Eq. (4.5) in the spacelike diamond < z, z < 1, see Section 3.4 



We will find it convenient to use the coordinates a, b vanishing at the center of the diamond: 

1,1, 

z = - + a + b, z — - + a — b. 

The sum rule functions F^a,i in this diamond: 

1. are smooth; 

2. are even in both a and b, independently: 

F dAt i(±a, ±b) = F d) A,i{a, b) ; 

3. vanish on its boundary: 

F dAjl {±l/2,b) = F dA;l {a,±l/2) = 0. 



(4.6) 



(4.7) 



Properties 1,2 are shown in Appendix [Dj Property 3 trivially follows from the definition of 
Fd,A,i, since both terms in the numerator contain factors zz{l — z)(l — z). 

A consequence of Property 3 is that the sum rule can never be satisfied with finitely many 
terms in the RHS. 



4.1 The sum rule in the free scalar theory 

To get an idea about what one can expect from the sum rule, we will demonstrate how it is 
satisfied in the free scalar theory. In this case d = 1, and only operators of twist A — / = 2 are 
present in the OPE 0x0 [24J,[25]. These are the operators 



Oajoc <j>d ltl ...d m <f>+... (A = / + 2, Z = 0,2,4,...)- (4 

be given I 

invariance, nonzero value of this correlator is consistent with Eq. (|3.6|) only if / is even. 



12 A formal proof of this fact can be given by considering the 3-point function (cf>(x)(j)(— x)O( M )(0)) . By x — * —x 
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The first term shown in (4.8) is traceless by 0's equation of motion, but it is not conserved. The 
extra bilinear in <p terms denoted by . . . make the operator conserved for I > (in accord with 



the unitarity bounds (3.3)), without disturbing the tracelessness. Their exact form can be found 
e.g. in [23] . 

In particular, there is of course the dimension 2 scalar 



O 



2,0 



1 

7T 



where the constant factor is needed for the proper normalization. At spin 2 we have the energy- 
momentum tensor: 



4 ,2 oc <fid^d u (, 



The operators with / > 2 are the conserved higher spin currents of the free scalar theory. 

The OPE coefficients of all these operators (or rather their squares) can be found by decom- 
posing the free scalar 4-point function into the corresponding conformal blocks, Eq. (4.4). We 
have [25], [21]: 

^ = 2' +1 ||j (l = 2n). (4.9) 



Using these coefficients, we show in Fig. [3] how the sum rule (4.5), summed over the first 
few terms, converges on the diagonal z = z of the spacelike diamond. Several facts are worth 
noticing. First, notice that the convergence is monotonic, i.e. all -F^a,/ entering the infinite series 
are positive. This feature is not limited to the free scalar case and remains true for a wide range 



of d, A, I; it could be used to limit the maximal size of allowed OPE coefficients (see footnote 16). 



Second, the convergence is uniform on any subinterval z G [e, 1 — e], e > 0, but not on the 



full interval [0, 1], because all the sum rule functions vanish at its ends, see Eq. (4.7). Finally, the 
convergence is fastest near the middle point z = 1/2, corresponding to the center a = b = of the 
spacelike diamond. Below, when we apply the sum rule to the general case d > 1, we will focus 
our attention on a neighborhood of this point. 



5 Main results 



In this section we will present a derivation of the bound (1.4), based on the sum rule (4.5). We 



assume that we are given a unitary CFT with a primary scalar operator of dimension d > 1. 
We consider the 4-point function ((jxfxfxp) and derive the sum rule (4.5), where the sum is over all 
primary operators appearing in the OPE 
about these operators, such as 



x 



We will use only the most general information 
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13 Thc energy-momentum tensor T„ v , which is a spin-2 primary of dimension 4, has to appear in the OPE, with 
a known coefficient |25j = ^d/(3y/cr) depending on the central charge ct of the theory. However, we are not 
making any assumptions about the central charge and will not take this constraint into account. It may be worth 
incorporating such a constraint in the future, since it could make the bound stronger. From the point of view of 
phenomenology, estimates of the electroweak S-parameter prefer models with small number of degree of freedom, 
hence small ct- 
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z=z 



Figure 3: The RHS of the sum rule in the free scalar theory, summed over I < 
0,2,4,8,16 (from below up) and plotted for < z = z < 1. The asymptotic ap- 
proach to 1 (dashed line) is evident. Notice the symmetry with respect to z = 1/2, a 



consequence of (4.6). 



1. only the operators satisfying the unitarity bounds (3.3) may appear; 

2. their spins I are even; 

3. all the coefficients pa,i are non-negative. 



We will prove the bound (1.4) by contradiction, 
operators of dimension A 



Namely, we will show that if only scalar 



> f(d) are allowed to appear in the OPE, the sum rule cannot be 
satisfied no matter what are the dimensions, spins, and OPE coefficients of all the other operators 
(as long as they satisfy the above assumptions 1,2,3). Thus such a CFT cannot exist! In the 
process of proving this, we will also derive the value of f(d). 



5.1 Why is the bound at all possible? 

Let us begin with a very simple example which should convince the reader that some sort of bound 
should be possible, at least for d sufficiently close to 1. 

The argument involves some numerical exploration of functions i*d,A,i entering the sum rule 



(4.5), easily done e.g. with Mathematica. These functions depend on two variables z, z, but for 
now it will be enough to explore the case < z = z < 1, which corresponds to the point X2 lying 
on the diagonal x\ — x% of the spacelike diamond in Fig. [2] We begin by making a series of plots 
of Fd t A,i for I = 2, 4 and for A satisfying the unitarity bound A > I + 2 appropriate for these spins 
(Fig. [i]). The scalar case I = will be considered below. We take d = 1 in these plots. 

What we see is that all these functions have a rather similar shape: they start off growing 
monotonically as z deviates from the symmetric point z = 1/2, and after a while decrease sharply 
as z — > 0, 1. These charecteristics become more pronounced as we increase I and/or A. We invite 
the reader to check that, for d = 1, these properties are in fact true for all F^a,i for even / > 2 and 
A > I + 2. By continuity, they are also true for d = 1 + e as long as e > is sufficiently small 
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14 



One can check that they are true up to d ~ 1.12. For larger d, F^ 4 2 (z = 1/2) becomes negative 
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F dA j (d=l,l=2) 



F dA i (d=l,l=4) 




0.25 0.5 0.75 



z=z 




0.25 0.5 0.75 



z=z 



Figure 4: The shape of -F rfj A,z for d 
unitarity bound. 



1, I = 2,4 and several values of A satisfying the 



Mathematically, we can express the fact that -F^a,; is downward convex near z = 1/2 as: 

^W>0 at 3 = 2=1/2, 
/ = 2,4,6..., A>/ + 2, 
1 < d < 1 + e. 



(5-1) 



Even before addressing the existence of the bound, let us now ask and answer the following 
more elementary question: could a CFT without any scalars in the OPE <p x <p exist? Eq. (5.1) 
immediately implies that the answer is NO, at least if d is sufficiently close to 1. 



The proof is by contradiction: in such a CFT, the sum rule (4.5) would have to be satisfied 



with only I > 2 terms present in the RHS. Applying the second derivative to the both sides of 

z 



(4.5) and evaluating at z 



1/2, the LHS is identically zero, while in the RHS, by (5.1), we 



have a sum of positive terms with positive coefficients. This is a clear contradiction, and thus 
such a CFT does not exist. 

To rephrase what we have just seen, the sum rule must contain some terms with negative 
F"(z = 1/2) to have a chance to be satisfied, and by (5.1) such terms can come only from I = 0. 
Thus, the next natural step is to check the shape of Fd,A,i for I — 0, which we plot for several 
A > 2 in Fig. [5] We see that the second derivative in question is negative at A = 2 (it better 
be since this corresponds to the free scalar theory which surely exists!). By continuity, it is also 
negative for A near 2. However, and this is crucial, it turns positive for A above certain critical 
dimension A c between 3 and 4. It is not difficult to check that in fact A c ~ 3.61 for d near 1. 

We arrive at our main conclusion: not only do some scalars have to be present in the OPE, but 
at least one of them should have A < A c ! Otherwise such a CFT will be ruled out by the same 
argument as a CFT without any scalars in the OPE. In other words, we have just established the 
bound A min < A c for d near 1. 

Admittedly, this first result is extremely crude: for instance, the obtained bound does not 
approach 2 as d — > 1. However, what is important is that it already contains the main idea of the 
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F dA i (d=l,l=0) 




0.25 0.5 0.75 1 
Figure 5: Same as Fig. |4j for I = 0. 



z=z 



method which will be developed and used with increasing refinement below. This idea is that we 
have to look for a differential operator which gives zero acting on the unit function in the LHS of 
the sum rule, but stays positive when applied to the functions F^a,i in the RHS. 

Now, some readers may find it unappealing that the method as we presented it above seems to 
be heavily dependent on the numerical evaluation of functions a i and their derivatives. Do we 



have an analytical proof establishing e.g. the properties (5.1)? - a purist of mathematical rigor 
may ask. 

Partly, the answer is yes, since the asymptotic behavior of Fo,a,i for large A and/or I is easily 
accessible to analytical means (see Appendix D). These asymptotics establish Eq. (5.1) in the 



corresponding limit. On the other hand, we do not have an analytic proof of Eq. (5.1) for finite 



values of A and I. Notice that such a proof must involve controlling hypergeometric functions 
near z = 1/2. No simple general expansions of hypergeo metrics exist near this point (apart from 
the one equivalent to summing up the full series around z — 0). Thus we doubt that a simple 
proof exists. 

Nevertheless, and we would like to stress this, the fact that we can establish Eq. (5.1) only 



numerically (with analytic control of the asymptotic limits) does not make it less mathematically 
true! The situation can be compared to proving the inequality e < ir. An aesthete may look 
for a fully analytical proof, but a practically minded person will just evaluate both constants by 
computer. As long as the numerical precision of the evaluation is high enough, the practical proof 
is no worse than the aesthete's (and much faster). 

To summarize, we should be content that numerical methods allow us to extract from gen- 
eral CFT properties (crossing, unitarity bounds, conformal block decomposition, . . . ) precious 
information about operator dimensions which would otherwise simply not be available. 



5.2 Geometry of the sum rule 

To proceed, it is helpful to develop geometric understanding of the sum rule. Given d and a 
spectrum {A, /} of O G <fi x 0, and allowing for arbitrary positive coefficients Pa,u the linear 
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A m ; n 5* A c 

sum rule violated 



A m j n A c 

critical case 



A m i n ^ A c 

sum rule satisfied 



Figure 6: The three geometric situations described in the text. The thick black line 
denotes the vector corresponding to the function F = 1. 



combinations in the RHS of (4.5) form, in the language of functional analysis, a convex cone C in 



the function space {F(a,b)}. For a fixed spectrum, the sum rule can be satisfied for some choice 
of the coefficients if and only if the unit function F(a, b) = 1 belongs to this cone. 

Obviously, when we expand the spectrum by allowing more operators to appear in the OPE, 
the cone gets wider. Let us consider a one-parameter family of spectra: 



E(A r 



{A,/ | A>A min (/ = 0), A > 1 + 2 (/ = 2,4,6...)} 



(5.2) 



Thus we include all scalars of dimension A > A E 
the unitarity bounds. 



u and all higher even spin primaries allowed by 



The crucial fact which makes the bound (1.4) possible is that in the limit A 



— > oo the 

convex cone generated by the above spectrum does not contain the function F = 1! In other 
words, CFTs without any scalars in the OPE 0x0 cannot exist, as we already demonstrated in 



Section [5TT] for d sufficiently close to 1. 

As we lower A min , the spectrum expands, and the cone gets wider. There exists a critical 
value A c such that for A min > A c the cone is not yet wide enough and the function F = 1 is still 
outside, while for A min < A c the F = 1 function is inside the cone. For A min = A c the function 
belongs to the cone boundary. This geometric picture is illustrated in Fig. |6j 

For A min > A c , the sum rule cannot be satisfied, and a CFT corresponding to the spectrum 
S(A min ) (or any smaller spectrum) cannot exist. By contradiction, the bound (1.4) with f(d) = 
A c must be true in any CFT. The problem thus reduces to determining A c . 

Any concrete calculation must introduce a coordinate parametrization of the above function 
space. We will parametrize the functions by an infinite vector |F( 2m ' 2n )| of even-order mixed 
derivatives at a = b = 0: 



p(2m,2ri) 



d 2m d 2n F{ 



a, b) 



=b=0 



(5.3) 



Notice that all the odd-order derivatives of the functions entering the sum rule vanish at this point 



due to the symmetry expressed by Eq. (4.6): 



p(2m+l,2n) _ p(2m,2n+l) _ p(2m+l,2n+l) _ g 

The choice of the a = b = point is suggested by this symmetry, and by the fact that it is near 
this point that the sum rule seems to converge the fastest, at least in the free scalar case, see 

Fig- HI 
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The derivatives (5.3) are relatively fast to evaluate numerically. Presumably, there is also no 
loss of generality in choosing these coordinates on the function space, since the functions entering 
the sum rule are analytic in the spacelike diamond. 

In terms of the introduced coordinates, the sum rule becomes a sequence of linear equations 
for the coefficients p^,i > 0. We have one inhomogeneous equation: 

i = E^S3, (5.4) 

and an infinite number of homogeneous ones: 

0=Ep^S, (5.5) 



We have to determine if, for a given A min , the above system has a solution or not. It turns 



mm > 2 which is of interest for us, 



out that in the range d > 1 and A 
RHS of the inhomogeneous equation 15. 41). In such a situation, if a nontrivial solution of the 
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all Fff] > in the 



homogeneous system (|5.5|) is found, a solution of the full system (|5.4|), (|5.5|) can be obtained by a 
simple rescaling. 



Thus for our purposes it is enough to focus on the existence of nontrivial solutions of the 

Geometrically, this means that we are studying the projection of 



homogeneous system (5.5) 
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the convex cone C on the F^>°) = plane. This projected cone, which is by itself a convex cone, 
may occupy a bigger or smaller portion of the _F(°> ) = plane, or perhaps all of it. Each of the 3 
cases shown in Fig. [6] can be characterized in terms of the shape of the projected cone, see Fig. [7) 

• A m i n > A c : This case is realized if the opening angle of the projected cone is 'less than 
7r', so that the homogeneous equations have only the trivial solution pa.z = 0. In a more 
technical language, the 'opening angle less than 7r' condition means that there exists a half- 
plane strictly containing the projected cone. If we write the boundary of this half-plane as 
A = 0, the linear functional A will be strictly positive on the projected cone (except at its 
tip). 

• A min = A c : This case is realized if the opening angle of the projected cone is equal to 7r in 
at least one direction. In other words, the boundary of the projected cone must contains a 
linear subspace passing through the origin. This subspace will be spanned by the projections 
of the vectors saturating the sum rule; the vectors from the "bulk" cannot appear in the sum 
rule with nonzero coefficients. 

• A m i n < A c : In this case the projected cone covers the whole plane. 



"That is, unless the scries in the RHS of (5.4) diverge. However, this situation does not occur in practice. 
16 Eq. (5.4| can become useful when studying other questions. E.g., it suggests that arbitrarily large OPE 



coefficients may not be consistent with the sum rule. It would be interesting to establish a model-independent 
theoretical upper bound on the OPE coefficients, which could be a rigorous CFT version of the NDA bounds in 
generic strongly coupled theories. For CFTs weakly coupled to the SM ('unparticle physics' [57]) experimental 
bounds on these OPE coefficients (also known as 'cubic unparticle couplings') have been recently considered in 
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Figure 7: The shape of the projected cones in each of the 3 alternative cases described 
in the text. 



Using this classification, we are reduced to studying the boundary of the projected cone. 
For practical reasons we will have to work with finitely many derivatives, i.e. with a finite- 
dimensional subspace of the function space or, equivalently, with a finite subset of the homogeneous 



system (5.5). The above geometric picture applies also within such a subspace. Satisfaction of 



the sum rule on a subspace gives (in general) weaker but necessary condition, so that we still get 
a valid bound (1.4) with f(d) = A c . As we expand the subspace by including more and more 
derivatives, the critical scalar dimension A c will go down, monotonically converging to the optimal 
value corresponding to the full system. 



5.3 Warmup example: d = 1 



Let us use this philosophy to examine what the sum rule says about the spectrum of operators 
appearing in the 0x0 OPE when has dimension d — 1. Of course we know that d — 1 
corresponds to the free scalar, see (1.2), and thus we know everything about this theory. In 
particular, we know that only twist 2 operators appear in the OPE, see Section 4.1 Our interest 



here is to derive this result directly from the sum rule. We expect the sum rule based approach 
to be robust: if we make it work for d = 1, chances are it will also give us a nontrivial result for 



d > 1. In contrast, the standard proof of (1.2) is not robust at all: it is based on the fact that 
the 2-point function of a d = 1 scalar is harmonic, and can hardly be generalized to extract any 
information at d > 1. 



The attentive reader will notice that our considerations from Section |5.1| were equivalent to 
retaining only the first equation out of the infinite system (5.5). This truncation did not control 
well the d 
A 



1 limit, since the obtained value of A r 



3.61 was well above the free theory value 
2. The next natural try is to truncate to the first two equations in (5.5). As we will see now, 



this truncation already contains enough information to recover the free theory operator dimensions 
from d — 1. 

Following the discussion in Section 5.2 we consider the projected cone — the cone generated 
by the vectors F = i*i,A,J projected into the two-dimensional plane (F^ 2 '°\ F^ 0,2 ^). For each 
I = 0, 2, 4, ... we get a curve in this plane, starting at the point corresponding to the lowest value 



of A allowed by the unitarity bound (3.3), see Fig. [8j It can be seen from this figure that: 
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Figure 8: The sum rule terms i*iA,z m the (_F( 2 '°) ; _F( > 2 )) plane. The shown curves 
correspond to / = 0, 2,4, 6. The arrows are in the direction of increasing A. The I = 
curve starts at A = 1 (A < 1.01 part is outside the plotted range); the I = 2,4,6 
curves — at A = / + 2. For large A the curves asymptote to the positive F^ 2 ' ^ axis, 
see Appendix |D} The shaded half-plane is the projected cone for a spectrum which 
includes the A = 2 scalar. 



1. the vectors corresponding to the twist 2 operators A = / + 2 lie on the line F^- 2 ' ^ = _F(°' 2 )p] 
while all the other vectors lie to the right of this line; 

2. the / = 0, A = 2 vector points in the direction opposite to the higher-spin twist 2 operators. 

The boundary of the projected cone is thus given by the line F^ 2 ' ^ = F^ ' 2 ^ if the spectrum 
includes the A = 2 scalar and at least one higher-spin twist 2 operator (e.g. the energy-momentum 
tensor). Otherwise the boundary will be formed by two rays forming an angle less than 7rJ^] By the 



classification of Section 5J2 it is only in the former case that the sum rule can have a solution. This 
case corresponds to A min = A c : the boundary of the projected cone contains a linear subspace 
passing through the origin. Thus we also have additional information: only the vectors from 
the boundary, i.e. those of the twist 2 operators, may be present in the sum rule with nonzero 
coefficients. 

The above argument appealed to the geometric intuition. For illustrative purposes we will also 



give a more formal proof. Taking the difference of the two first equations in (5.5), we get: 

o = E^( f S-^aO' ^ ,£0 - 

By property 1 above, for d — 1 all the terms in the RHS of this equation are strictly positive unless 
A = / + 2. Thus, only twist 2 operators may appear with nonzero coefficients. End of proof. 

It is interesting to note that in Fig. [8] the I = curve is tangent to the line F^ 2 ' ^ = F^ ' 2 ^ at 
A = 2.^ Were it not so, we would not be able to exclude the existence of solutions to the sum 
rule involving scalar operators of A < 2. 



17 This fact is easy to check analytically using the definition of Fd.A.i at d = 1. 

18 We ignore such subtleties as the possibility of a continuous scalar spectrum ending at A 
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Figure 9: The analogue of Fig. [8j for ci = 1.05. In this plot we started the I - 
at A = 2. The green line is the boundary of the projected cone for A min = A 
Fig " 



curve 
^ 3, see 

I0| The slope of this line is determined by the energy-momentum tensor vector. 



To conclude, we have shown that the spectrum of operators appearing in the sum rule, and 
hence in the OPE, of a d = 1 scalar consists solely of twist 2 fields and that, moreover, a A = 2 
scalar must be necessarily present in this spectrum. 

Isn't it amazing that we managed to find the whole spectrum using only the first two out of the 



infinitely many equations (5.5)? One may ask if by adding the complete information contained in 
the sum rule, a stronger result can be proved, namely that the full 4-point function of an arbitrary 
d — 1 scalar is given by the free scalar theory expression. This would constitute an independent 
proof of the fact that a d — 1 scalar is necessarily free. Such a proof can indeed be given [29], but 
we do not present it here since it is rather unrelated to our main line of reasoning. 

5.4 Simplest bound satisfying f(l) = 2 



We will now present the simplest bound of the form (1.4) which, unlike the bound discussed in 



Section 5.1 approaches 2 as d — > 1. The argument uses the projection on the (F^ 2 '°\ F^ 0,2 ^) plane 
similarly to the d = 1 case from the previous section. Since that method gave A min = 2 for d = 1, 
by continuity we expect that it should give A min ~ 2 for d sufficiently close to 1. 

To demonstrate how the procedure works, we pick a d close to 1, say d = 1.05, and produce 
the analogue of the plot in Fig. [HJ see Fig. [9| We see several changes with respect to Fig. [8j The 
energy-momentum tensor determines one part of the projected cone boundary (the green line), 
while the spins / = 4, 6, ... lie in the bulk of the cone. Continuation of the green line to the other 
side of the origin intersects the I = curve at the point corresponding to A = A c ~ 3. In the 



terminology of Section 5.2, this gives the critical value of A min . Namely, if A min > A c in the 



spectrum (5.2), the projected cone will have an angle less than n and the sum rule will have no 



solutions. On the other hand, for A min < A c the projected cone covers the full plane, see Fig. 10 



and a nontrivial solution to the first two equations of the system (5.5) will exist. For A min = A c 
the projected cone covers the half-plane shaded in Fig. [9j 

One can check that the same situation is realized for any d > 1. In particular, the slope of the 
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A<A C > '' 



A>A C 

Figure 10: The relative position of the I = vectors (red) with respect to the energy- 
momentum tensor vector (blue, pointing to upper right) determines the shape of the 
projected cone, see the text. If the cone contains the blue vector and both dashed red 
vectors, it covers the whole plane by their convex linear combinations. 



critical cone boundary, described by the linear equation 

F( 2 '°) -X(d)F^ = 0, 
is always determined by the energy-momentum tensor: 

i?(2,0) 



X(d) 



2,4.2- 



(5.6) 



(5.7) 



Once \(d) is fixed, the critical value of A m ; n is determined from the intersection of the line ( |5.6 ) 
with the I = curve: 

- X(d)F^ =0, F — F dAcfi . (5.8) 

The I = 0, A > A c points then lie strictly inside the half-plane F^ 2 ' ) - \(d)F i0 > 2 *> > 0. For 
A min > A c the cone angle is less than ir, and the sum rule has no solution. Thus the bound (1.4) 
must be valid with f(d) 
In Fig. 



A,. 
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we plot the corresponding value of f(d) found numerically from Eq. (5.7), (5.8) 



denoted /2(d) to reflect the order of derivatives used to derive this bound. As promised, the free 
field theory value A = 2 is approached continuously as d — > 1. 

The asymptotic behavior of /2(d) for d — > 1 can be determined by expanding the equations 
defining A c in power series in d — 1 and A c — 2. We find: 



f 2 (d) =2 + 1 Vd^l + 0{d-l), 
7 = [2(K + 1)/3] 1/2 ~ 2.929, 
K = (192 In 2 - 133)- 1 . 

This asymptotics provides a good approximation for d — 1 < 10~ 3 , see Fig. 
dependence in (5.9) can be traced to the fact that for d = 1 the I = curve was tangent to the 
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(5.9) 



The square root 



projected cone boundary at A = 2. The bound of Fig. 11 will be improved below by taking more 
derivatives into account, however the square root behavior will persist (albeit with a different 
coefficient). 
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1.00 1.02 1.04 1.06 1.08 1.10 1.12 



d=m 



Figure 11: /2(d) = A c as determined by solving Eq. (5.8). We plot it only for d rather 
close to 1 because in any case this bound will be significantly improved below. The 



dashed line shows the asymptotic behavior (5.9), which becomes a good approximation 
for d < 1.001. 



5.5 Improved bounds: general method 



As we already mentioned in Section 5.2, the bound will improve monotonically as we include 



more and more equations from the infinite system (5.5) in the analysis, i.e. increase the number 



of derivatives _F( 2m > 2n ) that we are controlling. We thus consider a finite basis £>, adding several 
higher-order derivatives to the F^ 2 ' ^ and F^ ' 2 ^ included in the previous section: 



B = { F^ m ' 2m h = { _p( 2 '°) ; _p(°> 2 ) ; . \. 



(5.10) 



According to the discussion in Section |5.2[ we have to study the boundary of the projected cone 
in the finite-dimensional space with coordinates (5.10). The logic in principle is the same as in the 
previous section. We will have a family of curves corresponding to I — 0,2,4. . . generating the 
projected cone. As we lower A min in the spectrum (5.2), the projected cone grows. For A min < A c 



it will cover the whole space. However, in this many-dimensional situation it is not feasible to 
look for A c by making plots similar to Fig. [9j We need a more formal approach. 



Such an approach uses the language of linear functionals, already encountered in Section |5.2 
A linear functional A on the finite-dimensional subspace with basis B is given by 



A(F) = ^A 2m , 2n F( 2m ' 2 ™) 



(5.11) 



where \2m,2n are some fixed numbers characterizing the functional. They generalize the single 



parameter X(d) from Section 5.4 to the present situation. 



Using linear functionals, the two non-critical cases of Fig. [7] can be distinguished as follows: 



ArT-iiv. ^> A r 



A m i n <C A c 



there IS a functional A such that 

HF dAjl ) > for all A, / G S(A min ) 

there is NO functional A such that 
A(F dA /) > for all A, I 6 S(A min ) 



(5.12) 
(5.13) 
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A numerical procedure which for any given A min finds such a positive A or shows that a non- 
negative A does not exists will be explained below. Assuming that we know how to do this, 
determination of A c becomes an easy task. First, we bracket A c from above and below by trying 
out a few values of A min and checking to which of the two above sets, A min > A c or A min < A c , they 
belong. Second, we apply the division-in-two algorithm, i.e. reduce the length of the bracketing 
interval by checking its middle point, etc. This achieves exponential precision after a finite number 
of steps. 

We will now explain the numerical procedure. Let us begin with the non-negative functional 



defined by Eq. (5.13), and comment later about the strictly positive case of Eq. (5.12). Eq. (5.13) 



can be viewed as a system of infinitely many linear inequalities for the coefficients \2m,2n- The 
infinitude is due to three reasons: 

• there are infinitely many spins I; 

• for each spin / the dimension A can be arbitrary large; 

• the dimension A varies continuously. 

To be numerically tractable, this system needs to be truncated to a finite system, removing 



each of the three infinities. We do it by imposing inequalities in (5.13) not for all A, / 6 S(A min ) 
but only for a 'trial set' such that 

• only finite number of spins I < Z max are included; 

• only dimensions up to a finite A = A max are included; 

• A is discretized. 

To ensure that we are not losing important information by truncating at / max and A max , we 
include into the trial set the vectors corresponding to the large I and large A asymptotics of the 



derivatives. The relevant asymptotics have the form (see Eq. (D.4) in Appendix D): 



-,(2m,2n) 



const 



(2m + l)(2n + l) 



(2\/2Z(l + a;)) 2m+1 (2\/2Z) 



2n+l 



A-/ 



X 



I 



> 



(5.14) 



where a constant const > is independent of m and n. This asymptotics is valid for I 
x <C I fixed. 



oo, 



Upon truncation to the trial set, Eq. (|5.13|) becomes a finite system of linear inequalities, a 

It is thus possible to determine if a solution 



particular case of the linear programming problem 
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exists (and find it if it does) using one of several existing efficient numerical algorithms (see [30J). 
In our work we used the classic Simplex Method as realized by the LinearProgramming function 
of Mathematica. 



If, using the linear programming, we find that (5.13) truncated to the trial set has no solution 



then a fortiori the full non-truncated system has no solution, i.e. non-negative A does not exist. 



19 A general linear programming problem consists in minimizing a linear function of several variables subject to 
a set of linear constraints (equalities and inequalities). Our problem is a particular case when all constraints are 
inequalities and the function to be minimized is absent (or, equivalently, it is constant). 
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MFdAj) ( d fixed > 1 fixed ) 




Figure 12: Solid curve: schematic typical dependence of the functional A(i^,A,z) on A. 
The functional A satisfies (5.15) for all A from a discrete trial set, which includes the 
dimensions A& and Afc +1 = A& + 8 A. Yet the functional may become slightly negative 
for Afc < A < Afc +1 . Dashed curve: same for a functional corresponding to a smaller 
5A. The violation of (5.13) at intermediate values has disappered. 



Thus we can safely claim that the considered A m i n brackets A c from below: A m i n < A c . The accu- 
racy of this bracketing will increase as we include more dimensions in the trial set (e.g. decreasing 
the discretization step). 



Let us now consider bracketing from above, which requires finding a functional satisfying (5.12 ) 



First, just as for (5.13), we truncate to a trial set. We'd like to use linear programming methods, 



however these methods do not work for strict inequalities. Thus we strengthen > in (|5.12|) to 
> £, where £ is a fixed small positive numbe 

A(F dAl )>e. (5.15) 



Then we can use the linear programming to find a solution of the truncated system. 

Unlike in the case of bracketing from below, to claim that indeed A min > A c , we have to 



check that the found A does not violate (5.12) for A, / not included in the trial set. This will not 
happen for I > Z max and A > A max if we take these parameters sufficiently large and include the 
asymptotics. The functional A may however become slightly negative at intermediate A and I, 
as a consequence of discretizing A. This is not unexpected, since the rays (A, I) that determine 
the boundary of the cone and thus the functional are determined with a fuzziness proportional to 
the discretization step 8 A. However, this violation will disappear as 8 A — > (for a fixed e), see 
Fig. [121 By choosing the discretization step smaller and smaller, we will be able to bracket A c 



from above with an arbitrary desired precision. 



20t 



In principle, Eq. (5.151 may become too constraining if all components of the vector Fd,A,i are O(e), which may 



happen in the large A limit. In this case one simply needs to rescale Fd t/ \,i by a constant factor. Each function 
Fd.A,i determines a ray in the finite-dimensional space, and such a rescaling does not change the content of the 



original Eq. (5.121. In practice, however, we never had to do such a rescaling. 
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Figure 13: Log-log plot of the maximal allowed anomalous dimension of <ft 2 , fe(d) — 2, 
versus the anomalous dimension of 0, d — 1. The dots correspond to the entries of 
Table 111 while the curve is the approximation (15.171). 



5.6 Best results to date 



In our numerical work we explored functionals (5.11) with the leading a-derivative up to F^ 6 ' ^ 



and with various choices of subleading derivatives. We now present our best results, which were 
obtained using the full list of derivatives with 2m + 2n < 6: 

| F ( 2m , 2 n) | (2m> 2n) = (6) 0); (4; 2)) (2; 4); (0; 6); (4) 0); (2) 2); (0) 4)) (2; Q)) (Q) 2) | (5 16) 



The bound f(d) = fe(d) corresponding to this choice is plotted in Figs. [T] and 13 Numerical 
values for several values of d are given in Table [l] For each d we give the bound f§(d) and the 
coefficients A 2m 2 n of the functional used to obtain this bound. These functionals were found using 
the linear programming method as described in the previous section. However, to check that our 
bound is true, one does not need to know how we found these functionals; it is enough to check 
that they indeed satisfy Eq. (5.12) with A min = fe(d). 



We have also computed the bound fe(d) for other values of d and found that it changes 
smoothly, interpolating between the points given in Table [TJ In the considered interval of d, the 
anomalous dimension fe(d) — 2 turns out to be well approximated (within ~ 2%) by the formula 



(see Fig. 13) 



f 6 (d) -2~ 1.79Vd- 1 + 2.9(d - 1) (1 < d < 1.35). 



(5.17) 



While our method would give a bound also for d > 1.35, we did not explore this range. The reason 
is that for d > 1.33 our bound exceeds A = 4 and starts getting not very interesting, taking into 
account the phenomenological motivations from Section [2] 

The coefficients \2m,2n i n Table [T] have been rounded with 6 significant digits; we have checked 
that the resulting slight violation of (5.12) is very smalj^] The anomalous dimensions f§(d) — 2 
in Table [l] approximate the optimal values, attainable using the subspace (5.16), from above to 



21 Thc rounded functionals violate (5.12 1 by about 10 4 for a few isolated values of A, I. This should be compared 
to the typical 0(1 -i- 1000) range of A(F) away from these points. 
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1.00 1.05 1.10 1.15 1.20 1.25 1.30 1.35 



Figure 14: From this plot one can get an idea how the bound monotonically improves 
as more and more derivatives are taken into acccount in the infinite system (5.5). The 



black (upper) curve is the simplest bound of Section 5.4, obtained by using only the 
second derivatives. The next two curves (green and red) correspond to the two sets 



(5.19). The lowest-lying blue curve is our best current bound obtained using the set 
(5~16). 



within 1% (due to the finite accuracy of the division-in-two algorithm used to bracket A c ); they 
have been rounded up with 3 significant digits. 

In obtaining these results, we used the trial set in the sense of the previous subsection with 
/ max = 10, A max = 20 (I = 0), A max = I + 10 (I > 2). We discretized A with step 5A = 0.01, 
decreased to 5A = 0.0025 around a few critical dimensions where the functional approaches zero, 



as in Fig. 12 for A& < A < A^+i. To take into account the asymptotic behavior (5.14), we have 



included into the trial set vectors with 



p(2m,2n) 



[(2m + l)(2n 
0, 



1)(1 



\2nl 



2m + In = 6 , 
2m + 2n < 6 , 



(5.18) 



obtained from (5.14) by rescaling and taking the I 
varying from x = to 10 with 5x = 0.01. We set e — 10 



oo limit. The parameter x in (5.18) was 
4 in rt5.15ft. 



We expect that including more derivatives in the list (5.16) should somewhat improve the 



bound, especially for d close to the upper end of the considered interval. We have observed a 
similar improvement trying out the functionals with the same leading a-derivative (6, 0) as in 



(5.16), but with smaller sets of subleading derivatives: 



Set 6a: (2m, 2n) = (6, 0), (4, 0), (2, 0), (0, 2); 

Set 6b: (2m, 2n) = (6, 0), (4, 0), (2, 0), (0, 2), (0, 4), (0, 6) 



(5.19) 



For illustration, we plot the corresponding bounds in Fig. [I4J including also the simplest second- 
derivative bound from Section 15. 4[ 
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d-1 


fa -2 


^6,0 


^4,2 


^2,4 


^0,6 


' X 4,U 


^2,2 


An a 
' X U,4 


^2,0 


^0,2 


10 -5 


0.00573 


1. 


-0.97747 


-1.06327 


-0.047622 


-116.282 


277.34 


49.2726 


3344.18 


-7170.94 


10 -4 


0.0185 


1. 


0. 


-0.00052291 


-0.999251 


-153.677 


48.2058 


126.869 


7344.28 


-8370.35 




0.0593 


1. 


-1.24135 


-0.845225 


-0.0234949 


-101.866 


276.743 


36.8807 


2626.56 


-6384.32 


0.01 


0.207 


1. 


-0.738985 


0.839453 


-1.00868 


-104.348 


14.8941 


88.7308 


4989.57 


-5288.51 


0.02 


0.31 


1. 


-1.04266 


1.22926 


-1.01672 


-87.2012 


0. 


73.9343 


4384.97 


-4356.11 


0.03 


0.397 


1. 


-0.669013 


0.980721 


-1.02839 


-100.265 


-8.59119 


82.6461 


5012.32 


-4772.52 


0.04 


0.476 


1. 


-1.14492 


1.51543 


-1.03388 


-77.4978 


-20.9209 


63.1812 


4200.66 


-3696.39 


0.05 


0.548 


1. 


6.31335 


-5.62554 


-1.0635 


-221.053 


-276.119 


338.147 


11062.1 


-4888.67 


0.07 


0.678 


1. 


8.20962 


-7.39456 


-1.07537 


-236.106 


-425.566 


439.223 


13130. 


-3371.37 


0.1 


0.849 


1. 


10.4578 


-9.55059 


-1.08398 


-255.95 


-620.004 


579.352 


16504.5 


-1662.47 


0.15 


1.11 


1. 


11.981 


-10.4519 


-1.14201 


-261.146 


-768.116 


649.227 


19160.9 


-46.5554 


0.2 


1.36 


1. 


12.7909 


-10.5811 


-1.20555 


-259.763 


-863.149 


676.812 


20924.6 


1544.21 


0.25 


1.61 


1. 


14.283 


-11.2729 


-1.28025 


-262.93 


-1027.33 


746.012 


24108.2 


3944.39 


0.3 


1.86 


1. 


18.0218 


-14.0589 


-1.36917 


-281.038 


-1467.79 


996.762 


32918.8 


9435.72 


0.35 


2.1 


1. 


24.6535 


-19.3331 


-1.49588 


-292.263 


-2367.92 


1493.39 


51357.2 


20484. 



Table 1: For several d we give numerical values of fa = f&(d) appearing in the bound 
(1.4) and the coefficients \2m,2n of a functional satisfying Eq. (5.13) with A min = / 6 , 
see the text. The functionals are normalized via A6,o = 1- 



6 Comparison to known results 

To further test our method, in this section we shall compare our bound to the operator dimensions 
in calculable CFTs. Several nontrivial tests are offered by exactly solvable CFTs in 2D. We will 
discuss these examples in subsection 1. On the other hand there are fewer calculable examples 
in 4D, and they all turn out to satisfy our bound in a somewhat trivial way. The point is that 
our bound on A lies abundantly above the line A = 2d, and none of the calculable models is 
significantly above this line. For instance in supersymmetric gauge theories the operators whose 
dimension is exactly calculable are chiral. In that case the relation d = 2r/3 holds, with r 
representing the -R-charge. Then, given a chiral scalar operator of dimension d, additivity of 
.R-charge implies that <fi 2 has dimension 2d, precisely on the A = 2d line, and thus our bound is 
trivially satisfied. 

In the case of large N theories Green's functions factorize at leading 1/N order, implying 
A = 2d + 0(1 /N). This relation does not provide a stringent test of our result unless d — 1 < 
0(1/N), which corresponds to an elementary, free, field at leading order in 1 /N. This situation can 
potentially be realized in variants of the Belavin-Migdal-Banks-Zaks (BZ) fixed point [3Tj , [32J . The 
simplest case of a non-Abelian gauge theory with matter consisting of charged fermions obviously 
does not provide a nontrivial check of our bound. This is because the gauge invariant operator 
of lowest dimension is ipip with dimension already close to 3. In order to have a chance to find a 
model that nearly saturates our bound we must necessarily add a scalar gauge singlet field <p to 
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the theory and look for a new BZ fixed point 22 Consider then a BZ model based on gauge group 
SU (N) with N-p fermionic flavors in the fundamental coupled to (p. The Lagrangian is 



C = -l^G^G 1 " + 0q + y(f>qq + A0 4 . (6.1) 



By dialing N-p/N = 11/2 — e with 1/N < e <C 1 the 1-loop /3-f unction is small and is compensated 
by the 2-loop contribution which comes with opposite sign {b > 0) 

9 d 1 2 7 o 2 iV 2 o 2 iV , . 

GS/i(^ 3 87T 2 S7H 

One can then easily check that (3(y) and /9(A) possess, already at 1-loop order, non-trivial zeroes 
satisfying y 2 ~ A and y 2 /87i 2 ~ e/(NN F ) ~ e/7V 2 . Notice also that, for such value of y 2 , its con- 
tribution to the 2-loop gauge /5-function is subleading in 1/7V, and therefore does not significantly 
affect the location of the zero of /3(g). The anomalous dimensions are given by 

y 2 N 2 

7^ = d- 1 = ci— = a x e, (6.3) 
7 02 = A-2 = 2 70 + c 2 — = 2a ie + a 2 -^. (6.4) 

87T Z iV z 

Again our bound is largely satisfied, just because 7^ arises at leading nontrivial order, at 1-loop. 
Notice that our bound in the small 7^ region is roughly 7^,2 < c^/7^, c ~ 1.79: in order to saturate 
it, 7^2 and 7^ should respectively arise at 1- and 2-loop order. This is never going to be the 
case if <fr has Yukawa couplings to fermions, but it could in principle be so if <fi is only coupled 
via quartic scalar couplings, as these do not lead to wave function renormalization at 1-loop. It 
is however easy to see that also this option does not help us to produce a nontrivial saturation 
of our bound. Indeed the fixed point condition necessarily implies that should enter at most 
linearly in scalar quartic couplings with charged scalars, otherwise the beta function /5(A) for the 
self-coupling would be strictly positive. Then even if a fixed point existed, with such limited, just 
linear, interaction 7^ and 7^2 would vanish at 1-loop. Saturation of our bound would then require 
7^2 = 2-loops and 7^ = 4-loops which seems unlikely to happen. Notice also that <fi 2 will mix 
with other invariant bilinears constructed with the charged scalars, so the "dimension of <p 2 " here 
means the lowest eigenvalue of an in principle complicated matrix of anomalous dimensions. 

One possible conclusion from the above discussion is that in order to saturate our bound, even 
at 7^ near 0, we necessarily need a theory at small N. In 4D we unfortunately have no other 
examples to play with. One obvious next try is to consider fixed points in 4 — e dimensions. Even 
if our bound strictly applies only to 4D theories, the comparison to fixed points in 4 — e is almost 
compulsory. The result, as we now show, is partly encouraging and partly frustrating. Consider 
the O(N) theory in 4 — e with Lagrangian 



22 See Section 5.1 of [T7] for a related discussion. 
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As was first studied in (23], this model has a fixed point at A(iV + 8)/487r 2 = e. There are 
two operators playing the role of 2 , the singlet Os = (fi a (fia and the symmetric traceless tensor 
Ot = 4>a4>b — (^/N)S a b(4 > c4 > c)- The computation of anomalous dimension gives [33] 

= Cfcee + ^ = [l - ~) + ^ - ^ 6 , (6.6) 

A 5 = A free + 75 = (2 - e) + ^±|e , (6.7) 
A T = A free + lT = (2 - e) + , (6.8) 

where in brackets we have indicated the free field scaling dimensions in D = 4 — e dimensions, 
cifree and Af ree , for and 2 respectively. In analogy, and "naive continuity" , with our study in 4D 
we should compare the anomalous dimension of the composite and elementary fields. Indeed the 
anomalous dimension of <fi arises only at two-loops so that we have 7^,5 oc v /t^ like in our bound! 
One always has jt < Is an d the most interesting relation is that between jj- and 7^ 

4 . . 

7t= - ^ = c n ^. (6.9) 

For iV > 3, < 1.79 consistent with our 4D bound. On the other hand Ci,c 2 > 1.79, above 
our bound. It is not clear what to make of this apparent contradiction, given that our bound 
surely applies only to 4D while here we are discussing a theory in 4 — e. On one side one would 
be tempted to argue that our bound smoothly extends to 4 — e, namely 

70 2 < c {t)\/l4> (6-10) 



with c(e) well behaved near e = 0. If our 4D result is correct then this cannot be the case. Instead 
it is possible that the relation between 7^2 and 7^ away from e = is more complicated than our 
result. Indeed we can view our 4D result as a bound on 7^,2/^7^ at ^7^ ^> e. The full result 
could be 

7^ < ^Ai^/e), (6.11) 

where A(x) is a function which interpolates between our coefficient c ~ 1.79 at x — 00, and a 
larger coefficient at x — 0, with a crossover around x ~ 1. For instance: 



5c 

A(x) = c+-„ — -, 5c>0, (6.12) 
x z + 1 



could do the job. 



6.1 Bounds in 2D CFT and comparison with exact results 

A wealth of information accumulated about exactly solvable CFTs in 2D [15] allows for a nontrivial 
check of our method. Much of our discussion in Sections 3,4,5 carries over to 2D with minimal, 
simplifying, changes. In 2D CFTs, we must make distinction between the global conformal group 
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SL(2,C), and the infinite-dimensional Virasoro algebra of local conformal trasnformations, of 
which SL(2, C) is a finite-dimensional subgroup Virasoro algebra plays crucial role in solving 
these theories exactly, but it has no analogue in 4D, and the results have to be expressed in 
SL(2, C) terms to allow for comparison. 

When we speak about primaries, descendants, conformal blocks in 2D theories, we must specify 
with respect to which group we define these concepts, Virasoro or SL(2, C). The former is standard 
in the 2D CFT literature, while it is the latter that is directly analogous to 4D situation 



2H 



Every Virasoro primary is a SL(2,C) primary, but the converse is not true. E.g. the stress 
tensor in any 2D CFT is a Virasoro descendant of the unit operator. To find SL(2,C) primaries, 
we need to decompose the sequence of all Virasoro descendants of each Virasoro primary (the so 
called Verma module) into irreducible SL(2, C) representations. While this is possible in principle, 
it may not be straightforward in practice. Nevertheless we know that SL(2,C) primaries have 
dimensions of the form 

A 5 l(2,c) = A Vir + n, n = or n > 2, 
where Av; r is a Virasoro primary dimension, and n is an integer. This is because the Virasoro 
operators which are not in SL(2, C) raise the dimension by at least 2 units. 
The unitarity bound for bosonic fields in 2D is 

A> I, 

where I = 0,1,2,... is the Lorentz spin. The SL(2,C) conformal blocks in 2D were found in 
|25j. |26pj in the same coordinates as before we have 

9A,i(u,v) = ^-[fA+i(z)f A -i(z) + (z~ z)\ . (6.13) 



Using the unitarity bound, the known conformal blocks, and the sum rule (4.5), which is valid in 
any dimension, we can try to answer the same question as in 4D. Namely, for a SL(2, C) scalar 
primary of dimension d, what is an upper bound on the dimension A m j n of the first scalar 



operator appearing in the OPE 0x0? I.e. we want a 2D analogue of Eq. (1.4). Since the free 
scalar is dimensionless in 2D, the region of interest is d — > <if ree = 0. 

In Fig. [15] we show such a bound on A min as a function of d obtained using the second derivatives 



F( 2 '°\ F^ '^ (thus this bound is analogous to the simple 4D bound from Section 5.4). The 
dependence looks approximately linear: 

f(d) = ff D \d) ~ 0.53 + Ad (2D, 2nd derivatives) . 

Unlike in 4D, this simplest bound does not approach the canonical value zero as d — > df rec = 0; 
we do not know if this has any deep meaning. Improvements of this bound using the method 



discussed in Section 5.5 are possible; see below 



We will now see how the bound of Fig. [15] checks with the known operator dimensions and 
OPEs in solvable unitary CFTs in 2D. Our first example is the free scalar in 2D. This CFT 
contains the so called vertex operator primaries given by an exponential of the fundamental field: 



V a = e 1 ^, [V a ]=a 2 . (6.14) 



SL(2, C) primaries arc sometimes called quasi-primaries in the 2D CFT literature. 
24 In contrast, explicit expressions for Virasoro conformal blocks are not known in general. 
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f 2 i2D \d) 




Figure 15: The solid (blue) line represents the simplest upper bound, in an arbitrary 
2D CFT, on the dimension A min of the first scalar in the OPE O x O of a dimension 
d scalar with itself. The dots show the position of the minimal model OPEs 0x0 and 
ip xip (see the text) in this plane. The dashed line corresponds to the free theory OPE 



(6.14). The bound is respected in all cases. 



The basic OPE of V a with itself has the form: 

V a xV a = V 2a . 

Thus we have d = a 2 , A = 4a 2 , which gives the dashed line in Fig. 15p* below the bound. 

A more interesting example involves the minimal model family of exactly solvable 2D CFT. 
The unitary minimal models (see [U],[31]) are numbered by an integer m = 3, 4, . . ., and describe 
the universality class of the multicritical Ginzburg-Landau model: 

£ ~ (<90) 2 + A0 2m ~ 2 . (6.15) 

For m = 3, the Ising model is in the same universality class. The central charge of the model, 

6 



1 



m(m — 1) ' 



monotonically approaches the free scalar value Cf ree = 1 as m — > oo. Intuitively, as m increases, the 
potential becomes more and more flat, allows more states near the origin (c grows), and disappears 
as m — > oo (free theory). 

Minimal models are called so because they have finitely many Virasoro primary fields (the 
number of SL(2,C) primaries is infinite). All Virasoro primaries are scalar fields O rs numbered 
by two integers l<s<r<m — 1, whose dimension is 

(r + m(r - s)) 2 - 1 

A r,s = Z— 7 — T • 6.16 

2m(m + 1) 



2oc 



Strictly speaking the bound in Fig. 15 was derived for real fields. However, we can apply it to the real parts 



which satisfy the OPE Re V a x Re V a ~ 1 + Re V 2a - 
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The Oi l is the unit operator (Ax,i = 0), while the field = 02,2 has the smallest dimension 
among all nontrivial operators: 

^ = A 2 , 2 = A lV (6-17) 
2m(m + 1) 



This field is identified with the Ginzburg-Landau field in (6.15). For m = 3 we have A 2j 2 = 1/8, 
which is the spin field dimension in the Ising model. 

It is convenient to extend the Virasoro primary fields to a larger range 1 < r < m — 1, 
1 < s < m, subject to the identification 



(r, s) <-> (m — r, m + 1 — s) 



(6.18) 



The fusion rules, which say which operators appear in the OPE O riSl x O r2S2 (but do not specify 
the coefficients) can now be written in a relatively compact form: 



O riSl X O r2S2 ~ 



(6.19) 



n 


- r 2 \ 




- r 2 \ 




- s 2 


+ 1,1*1 


- s 2 \ 



. min(ri + r 2 — 1, 2m — 1 — r± — r 2 ) 



For any m, the fusion rules respect a discrete Z 2 symmetry 
where ± 



(6.20) 



v — l) s_1 for m odd, (— l) r_1 for m even (this choice is dictated by consistency with 
(6.18)). This symmetry corresponds to the cf) — > —<fi symmetry of the Ginzburg-Landau model; in 
particular <p = 2 2 is odd under (6.20). 

We are interested in OPEs of the form OxO~l + + ... where both O and O have small 
dimensions. Two such interesting OPEs are 

l2 



0X0~l + 0^ + ... 

?/>x?/>~l + V ;2 + -- -, ^ = O h2 , cL 



1 

2 



2(m + 1 



(6.21) 
(6.22) 



Here <fi 2 and ip 2 are just notation for the lowest dimension operators appearing in the RHS. Note 
that for m = 3 we have ip = <fi via (6.18). Using the fusion rule (6.19) and the operator dimensions 

(6.23) 



(6.16) it is not difficult to make identification: 

<> 2 = Gig, A- = 1 



m = 3: 



m > 3: (j) 2 = 3j3 , A 02 



(Ising) 
4 



V" = Oi, 3 , A^ 2 = 2 - 



m(m + 1) ' 
4 



(6.24) 



m + 1 

In particular, we see that the ip x ?/> OPE does not contain 2 , which is precisely the reason why 
we are considering it 



20 



26 In general, (j) 2 does not appear in the OPE O ryS x O r<s for r = 1 or s = 1. The operators ^ has the lowest 
dimension among all these fields. 
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We are now ready for the check. Operator dimensions in both OPEs (6.21) and (6.22) are 
subject to the bound of Fig. 15 where we marked the corresponding points up to m = 10. We see 
that the bound is respected in all the cases, although the Ising model point lies remarkably close 
to the boundary. We have tried to improve the bound of Fig. [15] by including more derivatives in 



the functional, according to the general method described in Section 5.5 Although we have seen 



some improvement for lower valued of d, there was practically no improvement around the Ising 
spin dimension d — 1/8, so that also the improved bound was respected. One could wonder if the 
fact that the Ising model (almost) saturates the bound has any special significance. 

We have searched for other exactly solvable 2D CFTs which could provide checks of our 2D 
bound. E.g. some interesting OPEs can be extracted from the WZNW models. However, as far 
as we could see, none of them come close to saturating the bound. 

In conclusion, we would like to mention that some bounds for dimensions of operators appearing 
in the OPE of two primaries in 2D CFTs were derived in the past by Lewellen [37J and Christe and 
Ravanini [38J. Those bounds were however of a different nature than our bound (1.4). Roughly, 
the Lewellen-Christe-Ravanini (LCR) bounds say that IF a primary appears in the OPE, its 
dimension is not bigger than a certain bound. This is of course not the same as our result, which 
says that a certain primary MUST be present in the OPE, with the dimension not bigger than a 
certain bound. 

The methods of LCR are based on studying the monodromy of the conformal blocks near their 
singularities in the complex plane, as opposed to the more detailed information about the shape 
and size of conformal blocks at intermediate regular points used by us. They have to make a 
crucial assumption that only a finite number of singularity types exist, which means that the total 
number N of primaries appearing in the OPE, or at least the total number of such primaries having 
different A mod 1, is finite. This assumption is realized in Rational CFTs, but not in general. The 
LCR bounds become increasingly weak for large N and disappear in the limit N — ► oo. Thus it is 
doubtful that such methods could be useful in our problem, since we would like to be free of any 
assumptions about the spectrum of higher primaries. 



7 Comparison to phenomenology 



In this section we will comment on the precise relation of our main result (1.4) with the phe- 



nomenological discussion of Section [2j That discussion led to the constraints (2.21), (2.22) on the 
dimension d of the Higgs field operator H and on the dimension As of the first singlet in the 
H x W OPE, denoted W H . Are there any low-dimension non-singlets in this OPE? 

The standard considerations related to the p-parameter lead us to assume the "custodial" 
SO (4) = SU(2) L x SU{2) R as the global symmetry of the CFT 
Higgs field: 

hx + ih 2 
h 3 + ih 4 
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The real components h a of the 



H 



form a multiplet of primary scalars in the fundamental of SO (A). Their basic OPE will in general 
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Our conclusions would however not be affected by assuming just SU(2) x U(l). 



38 



A F =40 TeV 



A F =400 TeV 




Figure 16: Solid red: The region of (d, Ag) plane disfavored by phenomenological 
considerations of Section [2] In this region one or both constraints (2.21), (2.22) are 
not satisfied. We consider two cases: Ap = 40 TeV (left) and Ap = 400 TeV (right). 
Checkered blue: the region of the same plane which would be excluded by our bound 



(1.4) with / = / 6 if we knew that Ag < Ay. 



have the structure 



h a {x)h b (0) 



1 



\x 



2<l 



(5 ab l + C s \x\ As 5 ab (rfH)(0) + C T \x\ AT T {ab) (0) 



Cj\x\ 2 afjW(0) + ...) 



(7.1) 

Here we indicated the possibility for two symmetry structures in the even-spin sector: along with 
the 50(4) singlet WH, a scalar T ab transforming as a traceless symmetric tensor, i.e. the (1, 1) 
representation, will in general be present. E.g. in free theory T ab = h a h b — (l/4)5 ab (h c h c ) . On the 
other hand, in the odd-spin sector the first contribution will be associated with the conserved, 
dimension 3, SO {A) current j}? b \ which is in the antisymmetric tensor representation (1, 0) + (0, 1). 
Its coefficient Cj will be related to the normalization of the SO (4) Ward identity. 

Clearly (7.1) looks more complicated than (1.3). However, if we choose a = b = 1 (say) in 



(7.1), the current and other odd-spin fields drop out due to antisymmetry, and we get an OPE of 
the form (1.3). We just need to identify <fi = h\, <p 2 = H^H or 7n depending on which of these 



two scalars has smaller dimension. Eq. (1.4) then implies 



min(Ag,A T ) < f(d). 

We get a bound on Ag only assuming Ag < At- This is not satisfactory; in fact the reverse 
Ag > A^ seems likely (as it happens in the Wilson-Fisher fixed points in 4 — e dimensions, see 
Section [6]). However, we have reasons to believe (or hope) that our bound on min(Ag,Aj') is 
close to being saturated by some existing 4D CFT. One reason is that the truncation discussed 



in Section 5.6 shows clear signs of convergence. The other reason is that in 4 — e we do have 
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l<j> 2 ~ as dictated by our bound at small 7^. Taking this into consideration it makes sense to 



compare our bound f(d) to the phenomenological constraints on Ag. We have done that in Fig. 16 



assuming two different constraints on the scale Ap, a weaker one Ap > 40 TeV and a robust one 
Af > 400 TeV. According to the previous reasoning, even for the less likely situation As < A T , 
these plots indicate that there exists space for relaxing the flavor problem. The optimal value for 
the dimension of the Higgs field turns out to be between 1.2 and 1.3. Is it possible that a situation 
like this will be realized in nature? Hints of the answer to this question may soon come with the 
LHC. 

In the meantime, we believe that it should be possible to disentangle the contributions of WH 
and T in (7.1) and obtain a bound on A5 free of the assumption that As < At- As the above 
discussion clearly shows, such a bound cannot be found by considering only diagonal 0x0 OPEs 
where is a fixed h a component. One should try to use additional information contained in the 
nondiagonal OPEs with a 7^ b, something which we did not do in this paper. As it is apparent 
from (7.1), the global symmetry current and higher odd-spin operators will generally contribute 
to these nondiagonal OPEs. More work is needed to determine if the contributions of the odd-spin 
operators can be controlled in a model-independent way. 



8 Discussion and Outlook 

In this paper we have shown that prime principles of Conformal Field Theory, such as unitarity, 
OPE, and conformal block decomposition, imply the existence of an upper bound f(d) on the 
dimension A min of the first scalar operator 2 in the OPE of a scalar of a given dimension d. 

We developed a method which allows numerical determination of f(d) with arbitrary desired 
accuracy. The method is based on the sum rule, a function-space identity satisfied by the conformal 
block decomposition of the 4-point function (0000), which follows from the crossing symmetry 
constraints. In practical application of the method the sum rule is Taylor-expanded: replaced by 
finitely many equations for the derivatives. The bound f(d) improves monotonically as more and 
more derivatives are included; see Figs. [T, 10 for the best current bound obtained using derivatives 



up to the 6th order, and a sequence of weaker bounds obtained using fewer derivatives. 

We have checked that our bound is satisfied, by a large margin, in all weakly coupled 4D CFTs 
that we are able to construct. We have also derived an analogous bound in 2D and checked it 
against exact 2D CFT results. Again, the bound is satisfied, and in a less trivial way than in 4D, 
since the Ising model almost saturates it. 

Our results open up several interesting research directions: 

1. It should be relatively straightforward to improve our bounds, both in 4D and in 2D, us- 
ing our method but including more derivatives in the analysis. These improved bounds 
should monotonically converge to the optimal bound, corresponding to the infinite number 
of derivatives; we can already see signs of such convergence in Fig. [TUj 

2. One should search for more examples of CFTs which come close to saturating the bound, 
especially in 4D. 

3. The Dolan-Osborn closed-form expressions for conformal blocks are available only in even 
dimensions (up to D — 6). It is important to find expressions in 3D, comparable in simplicity 
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to (3.15 ) and (6.13 ). Then one could derive an analogous bound in 3D and confront it with the 
operator dimensions of known 3D CFTs, such as the O(N) universality classes. Although 
these theories are not exactly solvable, rather precise estimates for critical exponents and 
operator dimensions have been obtained using e-expansion, high-temperature expansion, and 
Monte-Carlo simulations [35]. For example, in 3D Ising model we have 7^ = 0.0183(4) and 
7^2 = 0.412(1) [SB]. 

4. It would be interesting to understand what is the appropriate extrapolation of our bound to 
4 — e dimension. This should explain why the comparison with Wilson-Fischer fixed points 
in Section [6] was not perfect. 

5. A very important but difficult problem is to find a genuine generalization of our bound to 
the situation when the CFT has a global symmetry. The case of 5*0(4) symmetry, readily 
generalized to SO(N), has been discussed in Section [7j 
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A Reality property of Euclidean 3-point functions 

In this appendix we would like to briefly discuss the reality property of 3-point functions, which 
was used at some point in our discussion. First of all, one can always choose a basis where each 
operator corresponds to an hermitian operator in the Minkowski space description. We work in 
such a basis. The reality properties of the Euclidean n-point functions for such operators are 
quickly deduced by analitic continuation of the Minkowskian correlators at space-like separation. 
Consider indeed the 3-point function 

G a ^ a ^(x 1 ,x 2 ,x 3 ) = (0^(x 1 )0^(x 2 )Or(x 3 )) , (A.l) 

where aij collectively denote the spin indices (we only consider bosons). When x\ 2 , x 23 and x 3 i 
are spacelike, the operators commute by causality, thus implying that G is real 

G^^( Xl ,x 2 ,x 3 )* = (0^(x 3 )0^(x 2 )0^( Xl )) = (0?{ Xl )0?(x 2 )0^{x 3 )} = G otl,a ' 2,a:i (x 1 ,x 2 ,x 3 ) 

(A.2) 

Continuation to the Euclidean then amounts to 

<-£ ^. %% £ % v $ j£ iyk — — 1^ 2 ; 3) j 

O a _> (-i) n «0%, (A.4) 
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where n a is the number of indices in {a}. Using the above rules, analytically continuing G ai,a2,a3 
from the spacelike patch we find the Euclidean functions 



-n a3 Q. 



Ql, 02,03 / 



-IX 







Eli x Eli 



-IX 







E2i x E2i 



-IX 







EZ) X E3 



(A.5) 



Now, by Lorentz and translation invariance G ai ' a2 ' a3 (xi, X2, x%) depends on the invariants x\^ 
while the tensor indices are covariantly reproduced by combinations of x^ with the invariant 
tensors t]^ and e pupa . It is now evident that if e pupcr does not appear in G ai,a2,a3 then the factor 
^n^+riaj+nag j n e q_ (A.5) is exactly compensated by a factor (— i^ n °'i +n °'2 +na 3 from the coordinate 
dependence of the tensor structure. In such a situation Q^' a2 ' as j s therefore real. On the other 



hand, contributions to that are proportional to one (equivalently an odd number) power of 
e iiupcr are p Ure imaginary. For the case that interests us in which two operators, say 0\ and O2 
have zero spin and O3 = O^ 1 "^ , the tensor structure of G pl : (x\, X2, £3) can only involve x^, 
and rf v and thus the corresponding euclidean function must be real. On the other hand the 
three point function for vector fields admits in conformal field theory a contribution proportional to 
e v.vpo which i s precisely proportional to the triangle anomaly diagram [39]. Therefore the euclidean 
3-point function for vector fields in CFT is in general complex. To conclude we notice that one 
can simply reproduce the result just discussed by formally assigning the following transformation 
property to the invariant tensors 



s. 



fJ,V 



-fiupcr 



fiupcr 



(A.6) 



B Closed-form expressions for conformal blocks 



The Dolan-Osborn result (3.15) is crucial for us, and we would like to say a few words about how 
it is derived, following [26 28 The main idea is that the conformal block is, in a certain sense, a 



spherical harmonic of the conformal group. In particular, it satisfies an eigenvalue equation 

T> X1 , X2 CB = -c A>/ CB , c Ali = 1(1 + 2) + A(A - 4), (B.l) 
where T> is a second-order partial differential operator acting on the coordinates x\p- This operator 



encodes the action of the quadratic Casimir operator of the conformal group 29 

1 



C 



D 



1 



(P.K. + KM = L A L A , 



2 — — p.f — 2 
on the state 0i(xi)0 2 (a ; 2)|O). The defining equation is 

C ■ 0i(xi)0 2 O 2 ) = [La, [L A ,(p 1 (x 1 )cf)2(x2)}} = V^^^x^fofa). 



On the other hand, ca,i in (B.l ) are nothing but the eigenvalues of the Casimir acting on conformal 
primaries [T2] : 

C-O M (0) = -c A , l O M (0). 



28 The first derivation [5S] is a brute force resummation of contributions of all conformal descendants of O and 
is not particularl y enl ightening. 

29 Here and in (B.2| we use the sign conventions of |14j for the conformal generators. 
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Since C is a Casimir, the same eigenvalue equation is simultaneously satisfied for all descendants 
of O. Thus Eq. (B.l) is a consequence of the OPE (3.9) and of the definition of conformal blocks. 

The explicit form of T> Xl>X2 can be found using the known expressions for the action of conformal 
generators on scalar primaries: 



[M Miy , (f>(x)] = i(xnd v - x v dn)<j> , 
[Kp, (f)(x)] = i(x 2 d^ - 2x^x ■ d - 2x l ^A (j 



(B.2) 



Eq. (B.l) can then be rewritten as a differential equation for g&j(u,v). The clever change of 
variables u,v — > z,z performed in (26] allows to find explicit solutions. In these variables the 
differential equation takes the form 



T> 1 



V z .- Z = z\\ - z)d 2 z + z\\ - z)d\ - (z 2 d z + z 2 d- z ) + 



2zz 



z — z 



[(l - z)d z - (i - z)d z ] 



The OPE fixes the asymptotic behavior: 

i M 



9a,i 



, A-i 
ZZ) 2 



2' 



z — z 



l+l 



z l+1 ), (z,z^0). 



With these boundary conditions the solution is unique and is given by (3.15). 



C z and z 



Here we comment upon the ranges of the variables z and z introduced in Section Since z and 
z are functions of the conformally-invariant cross-ratios, we can use conformal symmetry to fix 
some of the coordinate freedom. For instance, we can put three out of four points along a straight 
line, and moreover send one of them to infinity, as in Eq. (3.17). After that there still remains 



freedom to perform rotations leaving this line invariant, so that we can fix 

x 2 = (X 1 ,0,0,X i ). 
Now it becomes trivial to compute the cross-ratios: 



(C.1) 



u = x\ 2 = X\ + X 2 



r (.V, -l) 2 + Xl. 



Moreover, we can easily solve Eq. (3.16) for z and z: 



z = X\ + 1X4, z = z* . 



(C.2) 
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Thus we conclude that in the Euclidean, z is always the complex conjugate of z. Moreover, z is 
real if and only if all four points lie on a circle (which is a conformal image of the straight line in 
the above parametrization, see Fig. 17). 

Finally, Eq. (3. 18), (3. 19) are obtained from (C.1),(C.2) by Wick-rotating to the Minkowski 
time. 




Xo — X;\ 



x 1 



%3 



£4 



OO 



Figure 17: In Euclidean, z = z* . For the configuration of points chosen in this figure, 
z = X\ + iX^ is the complex coordinate of £2 in the X\ — X4 plane. 



D Asymptotic behavior 



In this Appendix we find large I and A asymptotics of derivatives of F d ^i at a = b = 0. It is 
useful to rewrite the definition of F d! A,i as follows: 



F d ,A,i(a,b) = h d (a,b) 
where we introduced the functions 



9dAA a M ~ 9dA,i{~ a M 



(D.1) 



S,w = [(i-z)(i-*)n/A,i 

h d (a,b) 



(zz) d - [(l-z)(l -z)] d ' 
These functions are smooth in the spacelike diamond. Moreover, it is not difficult to see that 

g(a, -b) = g(a, b), h d (±a, ±b) = h d (a, b) . 



In particular, from (D.l) we see the property (4.6) from Section 4 



Let us introduce the parameter 

8 = A - I - 2. 

As we will see below, there are three relevant asymptotic limits to consider: 
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• I large, 8 = 0(1) ; 

• I large, 8 large, 8 <C I 2 ; 

• 5 large, 8 ^> I 2 . 

In all these cases the large asymptotic behavior of derivatives will come from differentiating 
9a,i, which we write in the form 

9a,i = const(-) 1 — [k 2l+5+2 {z)k 6 (z) - (b -> -&)] . (D.2) 

In this Appendix by const we denote various positive constants which may depend on d, 8 or 
/ but are independent of the derivative order d 2m d 2n . These constant factors are irrelevant for 
controlling the positivity of the linear functionals defined on the cones. 

Starting from the following integral representation for the hypergeometric function (see [ID] ) 

T(r) f°° 

2 F 1 {a,b,c;x) = ) ' - \ e _w (l - e"*)^ 6- ^! - xe'^dt (Re c> Re b > 0) 



T(b)T(c-b)J 

and using the steepest descent method, we derive the large (3 asymptotics: 
k p {x) = e W2)h{x) [q(x) + 0(l/(3)}, 
M.) = lnf 4(1 -^n , g (x) 



x J ' 2 (l - V / T^r) <fl^x ' 

The leading asymptotic behavior appears when all the derivatives fall on the exponential factors 
in fp containing large exponents. Various prefactors appearing in (D.l) and (D.2) are not differ- 
entiated in the leading asymptotics. However, the a -1 and 6 _1 factors are responsible for changing 
the order of the needed derivative, as follows: 

( 2m ,2n) Const \(2m+l,2n) 

COnst(-Y (exp A) ( 2m+ l,2n + l) ; 



(2m + l)(2n + 1) 
A = (I + 8/2)h(l/2 + a + b) + (8/2)h(l/2 + a-b). 

To find the leading asymptotics, we expand A near a = b = 0: 

A — (1 + 8) [h(l/2) + a h'(l/2)] + lb h'{l/2) + {5/2)b 2 h"{l/2) + ..., (D.3) 
h'(l/2) = 2V2, h"{l/2) = -2V2. 



In the case 8 <C I the last term in (D.3) plays no role, and we get: 

,(2m,2n) COnst(-) 1 o 1 2m+l[.//i i \ n 2n+l 



F ^' 1 ~ ( 2 m + l)(2n + l) [/l/(1/2)(/ + 5 )] 2m+1 ^(V2)Z] 2n+1 , * < / a . (D.4) 
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This asymptotic is applicable for I large, while 5 can be small or large, as long as the condition 
5 ^ I 2 is satisfied; i.e. it covers the first two cases mentioned above. If on the other hand 5 ^> I 2 , 
it is the last term in (D.3) which determined the asymptotics of 6-derivatives, and we get 



FlXf n) ~ const(-y^^[h\l/2)5} 2 ^[h'\l/2)5r, 6 » I 2 . 
Because h"{l/2) < 0, the last asymptotics changes sign depending on the parity of n. 

References 

[1] G. Mack, "All Unitary Ray Representations Of The Conformal Group SU(2,2) With Positive 
Energy," Commun. Math. Phys. 55, 1 (1977). 

[2] M. A. Luty and T. Okui, "Conformal technicolor," JHEP 0609, 070 (2006) [arXiv:hep- 
ph/0409274]. 

[3] M. J. Strassler, "Non-supersymmetric theories with light scalar fields and large hierarchies," 
arXiv:hep-th/0309122. 

[4] S. L. Glashow, J. Iliopoulos and L. Maiani, Phys. Rev. D 2, 1285 (1970). 

[5] S. Dimopoulos and L. Susskind, Nucl. Phys. B 155, 237 (1979); E. Eichten and K. D. Lane, 
Phys. Lett. B 90, 125 (1980). 

[6] B. Holdom, Phys. Rev. D 24, 1441 (1981) 

[7] A. G. Cohen and H. Georgi, Nucl. Phys. B 314, 7 (1989). 

[8] L. Randall and R. Sundrum, Phys. Rev. Lett. 83 (1999) 3370 [hep-ph/9905221]. 

[9] W. D. Goldberger and M. B. Wise, Phys. Rev. Lett. 83 (1999) 4922 [hep-ph/9907447]. 

[10] J. Maldacena, unpublished; E. Witten, Talk at ITP conference New Dimensions in Field The- 
ory and String Theory, Santa Barbara http://www.itp.ucsb.edu/online/susy_c99/discussion; 
S. S. Gubser, "AdS/CFT and gravity," Phys. Rev. D 63, 084017 (2001) [hep-th/9912001]; 
H. Verlinde, "Holography and compactification," Nucl. Phys. B580, 264 (2000) [hep- 
th/9906182]. N. Arkani-Hamed, M. Porrati and L. Randall, "Holography and phenomenol- 
ogy," JHEP 0108, 017 (2001) [arXiv:hep-th/0012148]. R. Rattazzi and A. Zaffaroni, "Com- 
ments on the holographic picture of the Randall- Sundrum model," JHEP 0104, 021 (2001) 
[arXiv:hep-th/0012248]. 

[11] G. Mack and A. Salam, "Finite component field representations of the conformal group," 
Annals Phys. 53, 174 (1969). 

[12] S. Ferrara, R. Gatto and A. F. Grillo, "Conformal algebra in space-time and operator product 
expansion," Springer Tracts Mod. Phys. 67, 1 (1973). 

[13] I. T. Todorov, M. C. Mintchev and V. B. Petkova, "Conformal Invariance In Quantum Field 
Theory", Pisa, Italy: Sc. Norm. Sup. (1978) 273p 

[14] E. S. Fradkin and M. Y. Palchik, "Conformal quantum field theory in D-dimensions," Dor- 
drecht, Netherlands: Kluwer (1996) 461 p. (Mathematics and its applications. 376) 



46 



[15] P. Di Francesco, P. Mathieu and D. Senechal, "Conformal Field Theory," New York, USA: 
Springer (1997) 890 p 

[16] O. Aharony, S. S. Gubser, J. M. Maldacena, H. Ooguri and Y. Oz, "Large N field theories, 
string theory and gravity," Phys. Rept. 323, 183 (2000) [arXiv:hep-th/9905111]. 

[17] B. Grinstein, K. Intriligator and I. Z. Rothstein, "Comments on Unparticles," arXiv:0801.1140 
[hep-ph] . 

[18] J. Polchinski, "String theory. Vol. 1: An introduction to the bosonic string," Cambridge, UK: 
Univ. Pr. (1998) 402 p 

[19] M. Luscher, "Operator Product Expansions On The Vacuum In Conformal Quantum Field 
Theory In Two Space-Time Dimensions," Commun. Math. Phys. 50, 23 (1976). 

[20] G. Mack, "Convergence Of Operator Product Expansions On The Vacuum In Conformal 
Invariant Quantum Field Theory," Commun. Math. Phys. 53, 155 (1977). 

[21] K. G. Wilson and W. Zimmermann, "Operator Product Expansions And Composite Field 
Operators In The General Framework Of Quantum Field Theory," Commun. Math. Phys. 
24, 87 (1972). 

[22] G. Mack, "Duality In Quantum Field Theory," Nucl. Phys. B 118, 445 (1977). 

[23] A. Mikhailov, "Notes on higher spin symmetries," arXiv:hep-th/0201019. 

[24] B. Schroer and J. A. Swieca, "Conformal Transformations For Quantized Fields," Phys. Rev. 
D 10, 480 (1974). 

[25] F. A. Dolan and H. Osborn, "Conformal four point functions and the operator product 
expansion," Nucl. Phys. B 599, 459 (2001) [arXiv:hep-th/0011040]. 

[26] F. A. Dolan and H. Osborn, "Conformal partial waves and the operator product expansion," 
Nucl. Phys. B 678, 491 (2004) [arXiv:hep-th/0309180]. 

[27] H. Georgi, "Unparticle Physics," Phys. Rev. Lett. 98, 221601 (2007) [arXiv:hep-ph/0703260]. 

[28] J. L. Feng, A. Rajaraman and H. Tu, "Unparticle Self-Interactions and Their Collider Impli- 
cations," Phys. Rev. D 77, 075007 (2008) [arXiv:0801.1534 [hep-ph]]. 

[29] We are grateful to Gergely Harcos who kindly provided us with a key element of the mentioned 
proof. 

[30] W.H.Press, S.A.Teukolsky, W.T.Vetterling, B.P.Flannery, "Numerical Recipes. The Art of 
Scientific Computing", 3rd Edition (2007), Cambridge University Press. 

[31] A. A. Belavin and A. A. Migdal, "Calculation of anomalous dimensions in non-abelian gauge 
field theories," Pisma Zh. Eksp. Teor. Fiz. 19, 317 (1974); JETP Letters 19, 181 (1974). 

[32] T. Banks and A. Zaks, "On The Phase Structure Of Vector-Like Gauge Theories With Mass- 
less Fermions," Nucl. Phys. B 196, 189 (1982). 

[33] K. G. Wilson, "Quantum field theory models in less than four-dimensions," Phys. Rev. D 7, 
2911 (1973). 

[34] J. Polchinski, "String theory. Vol. 2: Superstring theory and beyond," Cambridge, UK: Univ. 
Pr. (1998) 531 p 



47 



[35] J. Zinn- Justin, "Quantum field theory and critical phenomena", Chapter 25. 

[36] M. Hasenbusch, K. Pinn, S. Vinti, "Critical Exponents of the 3D Ising Universality Class 
From Finite Size Scaling With Standard and Improved Actions", arXiv:hep-lat/9806012 

[37] D. C. Lewellen, "Constraints For Conformal Field Theories On The Plane: Reviving The 
Conformal Bootstrap," Nucl. Phys. B 320, 345 (1989). 

[38] P. Christe and F. Ravanini, "A New Tool In The Classification Of Rational Conformal Field 
Theories," Phys. Lett. B 217, 252 (1989). 

[39] E. J. Schreier, "Conformal symmetry and three-point functions," Phys. Rev. D 3, 980 (1971). 

[40] Bateman Manuscript Project, "Higher transcendental functions", vol. I, McGraw-Hill Book 
Company (1953). 



48 



